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APPENDIX  A.4 
AIRCRAFT 


A.4.1  SAMPLE  CALCULATIONS 

Section  4 presented  the  method  used  in  designing  the  nuclear  powered  aircraft.  A 
numerical  example  of  calculations  is  now  shown.  The  specifications  are  as  follows' 

a)  GW  = 2.000,000  lbs 

b)  W/S  = 60  Ibs/sq  ft. 

c)  t/c  = 0.18  at  the  wing  root 

d)  A = 0,4 

e)  AR  = 9 

f)  Runway  length  = 10,000  f*. 

g)  Fuselage  = 40  ft  x 428  ft 

h)  clmax  = 2 (FlaPs) 

i)  CLMAX  = 1 (No  flaps) 

j)  a = 0.926 

k)  LF  = 2.5 

l)  ra  = allowable  stress  for  2024  aluminum 

m)  Ec  = modulus  of  elasticity  in  compression 
Ec  = 10.7  x 106  psi  for  2024  aluminum 

Ec  = 29.0  x 106  psi  for  4130  steel 

n)  Mmac  = bending  moment  at  wing  MAC 

= 21 ,800,000  ft-lb 

o)  PMAC  = shear  load  at  wing  MAC  = 355,100  lbs 

p)  h1MAC=  height  of  beam  envelope  at  wing  MAC  = 8.06  ft 

q)  '1MAC  = chordwise  length  of  beam  envelope  at  wing  MAC  = 34.15  ft 

A.4.1 .1  STRUCTURAL  WEIGHT  (SECTION  4.1.1): 


By  Eq.  4.1. 1-1 


Structural  weight  = I 0.16  GW  + 


[■ 


9.8  GW 


] K. 


(100/W/S)0-63  W/S  J t/c 


K*KX 


Kt/C=  0.86 

= 0.95  (Ref.33,  Ch,  2,  p.  32) 

Kx  = 1.00 

therefore:  Structural  Weight  = 455,000  lbs 


1 


i 


L 


A.4-1 


A. 4. 1.2  TAKEOFF  THRUST  (SECTION  4.1.2.1): 

Takeoff  distance  = 10,000  ft 
From  Fig  4.1.2. 1-1,  K = 430 

Clto  = 0 75  Clmax  = 1,5 

Eq.  4. 1.2-1,  K = (W/S)(W/T)(1/CLT0)(1/«r) 

therefore:  W/T  = K CLT0  «r  /(W/S) 

W/T  = 9.9 

Available  W/T  = 8.91  (90%  of  W/T) 

Takeoff  thrust  = 225,000  lbs 

A.4.1.3  CHEMICAL  ENGINE  WEIGHT  (SECTION  4.1.2.2):  The  aircraft  must  have  the  capabil- 
ity to  complete  takeoff  if  one  chemical  engine  is  lost  during  takeoff  (Ref.  116,  p.  13). 

Use  six  50,000  lb  thrust  engines. 

Eq.  4.1. 2-2  Engine  weight  = EXP  [ I:.36  ] 

Engine  weight  --  8,333  Ibs/engine 

Total  weight  = 50,000  lbs 

A.4.1.4  CHEMICAL  FUEL  WEIGHT  (SECTION  4.1.3):  Fuel  weight  = 0.16  GW 
Fuel  weight  = 320,000  lbs 

A.4.1.S  FIXED  EQUIPMENT  WEIGHT  (SECTION  4.1.4):  Fixed  equipment  weight  = 0.10  GW. 
Fixed  equipment  weight  = 200,000  lbs 

A.4.1.6  L/D  CALCULATION  (SECTION  4.1.5.1):  Let  aircraft  speed  = 350  kts  (mission  require- 
ment) at  30,000  ft. 

p = 0.00089  slugs/cu  ft 

V = 591  ft/sec 


Sw  = GW/(W/S)  = 33,333  sq  ft 


* 


■i 

I 

i 


!• 


1 


n* 

' , 


( 

1 ! 


Eq  4. 1.5-3  CL  - 2 GW/(/.  Sw  V2) 

CL  « 0.386 

Eq.  4. 1.5-6  F * 129  4 0.007  Sw  + (0.0021  )(Number  of  engines)(T/engine)0,7 
From  A. 4. 1.3.  Six  50,000  lb  thrust  engines  required. 

F = 387. 

Eq  4. 1.5-5  C0p  = F/Sw 
CDp  = 0.0116 

Eq.  4.1. 5-7  CD|  - CLJ/(rr  AR  e) 

CD|  * 0.00659 

CDC  = 0.002  (Ref.  33,  Ch.  2,  p.  44) 

Eq.  4. 1.5-4  CD  = CDp  + CQ|  + CQC 
CD  = 0.02019 
Eq.  4.1. 5-2  L/D  = CL/CD 
therefore:  L/D  = 19.12 

A.4.1.7  NUCLEAR  REACTOR  SYSTEM  WEIGHT  (SECTION  4.1.5): 

By  Eq.  4. 1.5-1:  T = GW/(L/D)  and  Power  = TV 
From  A.4.1.6:  V = 581  ft/sec  and  UD  = 19.12 
therefore:  T = 105,000  lbs 

Power  =•  61,820,000  ft-lb/sec  = 84  mw 
At  35%  efficiency  and  15%  for  additional  power  (Section  4. 1.5.2) 

Reactor  power  = (Power)(1/.35)(1.15) 

Reactor  power  =-275.4  mw 

From  Figure  4.1 .5.2-1  (Gas  Reactor):  Reactor  system  weight  = 603,000  lbs 


^ - .. 
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A.4.1.8  PAYLOAD  (SECTION  4.1.6): 


Payload  = GW  - (Structural  + Engine  + Fuel  + Fixed  + Reactor) 


theref  jre:  Payload  = 372,000  lbs 


A.4.1.9  TAKEOFF  PERFORMANCE  (SECTION  4.4.1): 


By  Eq.  4.4.1-1  V «=  (2  GW/(p  ST  CLT0))V2 


Let  p = 0.00231  Slugs/cu  ft  (1,000  ft  field  elevation) 


clmax  ~ 2:  clto  ” 1,5 


V = 186  ft/sec  = 110  kts 


A. 4.1. 10  CLIMB  PERFORMANCE  (SECTION  4.4.2): 


R/C  = 600  fpm  = 10  ft/sec  (mission  requirement)  at  sea  level 


p = 0.002378  slugs/cu  ft;  CL  = 0.4  (angle  of  attack  function) 


By  Eq.  4.4.2-2:  TQ  = 


GW  R/C  (p  ST  CL)0  5 


(2  GW)°- 


therefore:  Tc  = 56,300  lbs  additional  thrust  required  at  sea  level 


At  service  ceiling  of  40,000  ft:  R/C  = 100  fpm 


p = 0.0005857  slugs/cu  ft 


By  Eq.  4.4.2-2:  Tr  = 4700  lbs  additional  thrust  required  at  service  ceiling 


A.4.1.11  CRUISE  PERFORMANCE  (SECTION  4.4.4): 


Maximum  L/D  = 20.30  (determined  by  iteration  of  L/D  equation  over  a range  of  V) 


Cruise  L/D  = 0.8  Maximum  L/D  = 16.24  (Refs.  9 and  94) 


By  Eq.  4. 1.5-1:  Cruise  thrust  = GW/Cruise  L/D 


Cruise  thrust  = 123,000  lbs 


By  Eq.  4.4.4-1 


)n'"- 

p ST  L/D  (CDp  + Cpg) 
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p = 0.000890  slugs/cu  ft  (30,000  ft) 


I 

I 


i 

I 


\ 


t 


therefore:  with  data  from  A.4.1.6:  V = 640  ft/sec  = 379  kts 


A.4. 1.11.1  STALL  SPEED/THRUST  (SECTIONS  4.4.4.1  AND  4.4.4.2): 


By  Eq.  4.4.1-1:  V = ( 


2 GW 

P ST  CLMAX 


p = .000890  slugs/cu  ft 


VSTALL  = 367  ,,/seC  = 217  kts 


By  Eq.  4.4.4-2:  Stall  UD 


^LMAX 


CDP  + CDC 


+ 


C 2 

°LMAX 
it  AR  e 


Stall  thrust  = GW/(Stall  UD)  = 115,600  lbs 

since  required  cruise  thrust  is  greater  than  required  stall  thrust. 

By  Eq.  4.4.2-3:  Total  thrust  = Climb  thrust  + MAX  (Cruise  thrust  or  Stall  thrust) 
therefore:  Total  thrust  = 138,000  lbs  at  30,000  ft 


A.4.1.11.2  MAXIMUM  SPEED  (SECTION  4.4.4.3): 

UD  at  VMAX  = GW/Tt  = 14.49 
therefore  by  Eq.  4.4.4-1  (above) 

V = 777  ft/sec  = 460  kts  = Mach  0.78 
A.4.1.12  DESCENT  PERFORMANCE  (SECTION  4.4.5): 
Maximum  R/D  = 8000  fpm  (mission  requirement) 
By  Eq.  4.4. 5-2 

„ 2 GW  ( R/D  ^ \ 

s P V2  CDS  sin  tf>  ' V L/D  7 

Cos  = 1 (Ref.  73;  Ch.  13,  p.  11) 

V = 640  ft/sec,  L/D  = 16.24  (From  A.4.1.11) 
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IP 


</»  = 60'  (Figure  4.4. 5-1) 

/>  = 0.000890  (30,000  ft) 

therefore:  Ss  = 3321  sq  ft  of  speed  brf-.Kes 

A.4.1.13  LANDING  PERFORMANCE  (SECTION  4.4.6): 

By  Eq.  4.4. 6-1: 

(V  2 _ v 2 \ 

“^Ig — ~ + 50  ) 

By  Section  A.4.1.6  with  CLMAX  = 2:  L/D  = 10.50 

By  Eq.  4.4.M:  VSQ  = [ 2 GW l(p  ST  CLMAX)  j V* 

Vso  = 166  ft/sec  = 98  kts  at  3000  ft  field  elevation 
Vso  = 1.3  V;  VTQ  = 1.15  V 

Vso  = 216  ft/sec  = 130  kts:  Vtd  = 191  ft/sec  = 113  kts 
therefore:  dso  = 2184  ft 
By  Eq.  4.4.6-2: 

Ground  deceleration  distance,  dG  = VTD2/2a 

a = 6 ft/sec2  (Ref.  33,  Ch.  2,  p.  14) 
dQ  = 4378  ft 

By  Eq.  4.4.6-3:  Landing  field  required  = (d50  + dG)1.67 
Landing  field  required  = 10,960  ft 
By  Eq.  4.4.6-4:  R/D  at  Touch  down  = VTD(50/dso) 
therefore  R/D  at  Touch  down  = 4.4  ft/sec  = 262  fpm 
A. 4.1. 14  AIRCRAFT  TIRES  (SECTION  4.3.7.1): 

Minimum  footprint  = GW/250  psi  = 8000  sq  in 


J 


i 
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With  56  x 16  tires  with  a 241  sq  in  footprint/tire 

Number  of  tires  - minimum  footprint  = 33  2 
footprint/tire 

therefore  use  36  tires 

therefore:  aircraft  footprint  = 8700  sq  in 

and  aircraft  pressure  print  = 230  psi 

A. 4.1. 15  LANDING  GEAR  STRUT  (SECTION  4.3.7.2): 

By  Eq.  4.3. 7*1 : Energy  = V2  m V2 

By  Eq.  4.3.7-2:  m = (GW/g)  = 62,112  lb  sec2/ft 

V = Max  R/D  at  Touch  down  = 500  fpm  (Aircraft  limitation) 

therefore:  Energy  = 2,156,600  ft-lbs 

with  four  main  gear:  Energy/Gear  = 539,165  ft*lbs 

with  Aluminum  7075  struts  with  3 in.  thick  walls  the  allowable  pressure  = 

1,800,000  Ibs/sq  ft  = 12,500  psi 

By  Eq.  4.3. 7-3:  Energy  absorbed  = Strut  pressure  change  x strut  volume  change 
or:  Energy/Allowable  pressure  = Volume  change 

therefore:  Volume  change  = 0.30  cu  ft  with  a strut  deflection  of  2 ft,  internal  strut  area 

is  0.15  sq  ft 

with  a triple  oleo  strut: 

internal  oleo  diameter  = 5.25  in. 

outside  oleo  diameter  = 11.25  in. 

A.4.1.16  AIRCRAFT  COMPONENT  SIZE  AND  WEIGHT:  An  initial  GC  was  located  at  260  ft 
from  the  canard  aircraft  nose  with  the  use  of  structural  component  weight  as  a percent  of  total 
structural  weight  as  shown  in  Table  A.4.1 
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TABLE  A.4.1  COMPONENT  WEIGHTS  (REF.  33,  CH.4,  P.  41) 


COMPONENT 

% OF  STRUCTURE 
WEIGHT 

LOCATION  OF  CG 

WING 

40.0 

25%  M.  A.  C. 

FUSELAGE 

30.6 

40%  OF  FUSELAGE 

TAIL 

7.8 

25%  TAIL  M.  A.  C. 

NACELLE 

5.6 

40%  OF  NACELLE 

LANDING  GEAR 

16.0 

AT  AIRCRAFT  CG 

Wing/Canard  area: 

By  Eq.  4.3. 1-1 : Sc  = 


Cj^yy  dW 

CLCdC~ 


where:  CLW  « wing  lift  coef. 


CLC  = canard  lift  coef.  = CLW 


dW  = distance  from  wing  quarter  chord  to  CG 
dC  = distance  from  canard  quarter  chord  to  CG 


therefore:  Sc  = 6683  sq  ft 

Both  wing  and  canard  have  lift  and  total  lift  area  required  is  33,333  sq  ft.  Total  available  lift 


area  is  Sw  + Sc. 


therefore:  reduce  Sw  to:  ( 


■'w 


®w  + ®c 


and  Sc  to: 


;cs 

+ s, 


) 


therefore:  Sw  = 27,800  sq  ft;  Sc  = 5600  sq  ft 
Vertical  tail  area: 


By  Eq.  4.3.1-2  SVT  = 


0.5  DF  CDF  (dF2  - dA2) 

CDVT  dVT 
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where:  SVT  = vertical  tail  area 

DF  = fuselage  diameter  = 40  ft 

CDF  = fuselage  drag  coefficient  = 0.5 

dF  = distance  fwd  of  CG 

dA  = distance  aft  of  CG 

Cdvt  = vertical  tail  drag  coef.  = 0.91 

dVT  = distance  from  CG  to  vertical  tail  quarter  cord 

therefore:  SVT  = 4,300  sq  ft 

Wing/canard  planform: 

By  Eq.  4.3.4-2:  Span,  bw=  (>H  Sw)  »/2 
By  Eq.  4 3.4-3 

Hoot  chord,  C„  - ( J^r x ) i^-)1 


By  Eq.  4.3  4-4:  Tip  chord,  CT  = X CR 
By  Eq.  4.3.4-5: 

MAC  = 2/3  CR  ( l-tX  t-i2  ) 
R ' 1 + X ' 


therefore; 


Wing 


Canard 


bw  (ft) 
CR  (ft) 
CT  (ft) 
MAC  (ft) 


500 

*225 

79 

35.6 

31.6 

14.2 

58.7 

26.5 

Component  weights:  (Section  4.3.2) 


By  Eq.  4.3.2-6  Wing  weight  = 1624  ( W/S  Sw  LF  bw  Sw 

V wing  root  thickness  x 109 


j 0.584 


Wing  weight  = 210,200  lbs 
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By  Eq.  4.3.2-7:  Canard  weight  = 1624  / W/S  Sc  LF  bc  Sc 


0.584 


canard  root  thickness  x 109 


LI 

it 


Canard  weight  = 31,200  lbs 


By  Eq.  4.3. 2-3  Vertical  tail  weight  = 1.197  (SVT)' 


24 


Vertical  tail  weight  = 38,000  lbs 


By  Eq.  4.3. 2-4:  Fuselage  weight  = 


°'273  t ( "height"”  ^ ° ^ ^maxJ0’,#  (wetted  area)1-05  ] 


0.88 


where:  QMAX  = Vt  i>  V2  and  V = 524  ft/sec  at  sea  level 


qmax  = 326  ,b/sc|  ft 


thus:  Fuselage  weight  = 188,000  lbs 


By  Eq.  4. 3.2-5  Landing  gear  weight  = 0.00916  (GW)1,124 


i 

* 

I • 

! 

i 

‘ 


Landing  gear  weight  = 111,000  lbs 

A.4.1.17  AIRCRAFT  WEIGHT  AND  BALANCE  (SECTION  4.3.8.1):  No  chemical  fuel  will  be 
used  in  this  example.  The  fuel  may  be  used  to  shift  the  CG  position  if  necessary  The 
components  weight  moments  are  summed  about  the  aircraft  nose  so  that: 


CG  moment  arm  - Sum  of  (component  weights  x moment  arms) 

total  weight 


A. 4-10 


therefore- 


ITEM 

WEIGHT  (lbs) 

LOCATION  (ft) 

MOMENT 
ft-lbs  x 106 

Wing 

210,200 

300 

63.06 

Canard 

31,200 

60 

1.87 

Vert.  Tail 

38,000 

360 

13.68 

Fuselage 

188,000 

172 

32.34 

Landing  Gear 

111,000 

300 

33.30 

Chemical  Engines 

50,000 

300 

15.00 

Fixed  Equipment 

200,000 

75 

15.00 

Nuclear  Reactor 

603,000 

300 

180.90 

Payload 

372,000 

200 

74.40 

Total 

1,803,400 

429.21 

Thus:  CG  at  238  ft  from  aircraft  nose  with  no  chemical  fuel.  Up  to  1 96,600  lbs  of  chemical 
fuel  may  be  carried  by  the  aircraft.  This  results  in  a tradeoff  between  the  payload 
and  chemical  fuel  weights  and  iteration  of  the  aircraft  design. 


A.4.1.18  WING  SECTION  PROPERTIES  (SECTION  4.3.4.5):  Section  property  calculations  are 
made  at  the  wing  station  for  the  mean  aerodynamic  chord.  The  beam  section  area  moment  of 
inertia  is  calculated  from  Eq.  4.3.4-8. 


(LF)Mhi  (2.5)  (21,800,000)  (8.06) 
2 Tg  ” (2)  (40,000)  (144) 


= 38.12  ft4 
From  Eq.  4.3.4-9 

hi-Vh'"17f  '</»» 

= 7.99  ft 

The  beam  cap  thickness  is  calculated  from  Eq.  4.3.4-10 


tc  = Vs  (hi  - h2)  = 8 062'  7 99  = 0.035  ft  = 0.42  in. 
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The  beam  shear  web  thickness  is  calculated  from  Eq.  4.3.4-12 


(LF)  P (2.5)  (355,100) 
w 2 h2  ra  (2)  (7.99)  (40,000)  (144) 


= 0.0096  ft  = 0.115  in. 


So  set  tw  = 0.125  in. 


A.4.1.19  RIB  SPACING  CALCULATION  (SECTION  4.3.4-6):  The  stiffener  spacing,  b.  is  found 
from  Eq.  4.3.4-13 

h _ ( KEc  V ) 1/2  [ (6.35)  (10.7  X 10")  (0.42)2  I1/2  „ „ 

“ L 40000 J -1715in- 


The  rib  spacing,  a,  is  equal  to  the  a/b  ratio  times  b. 
a = 4b  = (4)  (17.15)  = 68.6  in. 

The  number  of  ribs  in  the  wing  is  found  by  dividing  the  wing  span  by  the  rib  spacing. 


# ribs  = 


499  9 X 12 
68.6 


= 87.4  = 87 


A.4.1.20  DIVERGENCE  SPEED:  Divergence  speed  is  found  from  Eq.  4.3.4-14.  The  matrices 
[ejl/i-dyl.  [e]  and  | W|  are  presented  along  with  their  matrix  product.  Elements  of  each  matrix 
are  calculated  at  the  six  wing  stations  given  as  distances  in  ft  from  the  aircraft  centerline: 


yi  = 

241.43 

y*  = 

216.46 

ya  = 

176.74 

y < = 

124.98 

ys  = 

64.69 

y«  • 

0 

(cj  = 


l/J-  dyl  - 


33.36 

0 0 

0 0 0 

ft 

0 3812  0 

0 0 0 

0 

0 45.68  0 0 0 

0 

0 0 

55.54  0 0 

0 

0 0 

0 67.03  0 

0 

0 0 

0 0 79.35  J 

11.52 

11.52 

11.52  11.52 

11.52 

0 

11.52 

6.35 

6.35  6.35 

6.35 

0 

11.52 

6.35 

2.36  2.36 

2.36 

0 

11.52 

6.35 

2.36  .85 

.85 

0 

11.52 

6.35 

2.36  .85 

.26 

0 

0 

0 

0 0 

0 

0 

A.4-12 


[c]|/-J-dyl[e]|W|  = 


5.18 

0 

0 0 

0 

0 

0 

5.48 

0 0 

0 

0 

0 

0 5.96  0 

0 

0 

0 

0 

0 6.59  0 

0 

0 

0 

0 0 

7.31 

0 

0 

0 

0 0 

0 

8.09 

16.94 

0 

0 

0 

0 

0 

32.72 

0 

0 

0 

0 

0 

46.27 

0 

0 

0 

0 

0 

56.67 

0 

0 

0 

0 

0 

63.21 

0 

0 

0 

0 

0 

.0337 

.0689 

.1060 

.1435 

.1777 

.0385 

.0434 

.0667 

.0903 

.1118 

.0462 

.0520 

.0297 

.0403 

.0498 

.0562 

.0632 

.0362 

.0175 

.0217 

.0678 

.0763 

.0436 

.0212 

.0081 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

32.72 


The  eigenvalue,  r,  of  the  matrix  product  is  0 2934  x 1 0®.  The  divergence  speed  is  found  from  Eq. 
4.3.4-17. 


' 2 G -v  / (2)  (4  X 106)  (144) 

/»7Cl  v V (0.002378)  (0.2934  X 106)  (4.94) 


= 578  ft/sec  at  sea  level 
= 343  kts 

A. 4. 1.21  REACTOR  MOUNTING  STRUCTURE:  (SECTION  4.3.6):  The  design  algorithm  for 
the  largest  column  is  given  here.  It  has  an  axial  g-loading  of  10.61,  so  the  compressive  load  is 
(10.61)  (800,000)  = 8.488  x 10®  lbs  and  is  150  in.  long.  From  equation  4.3.6-1  the  cross 
sectional  area  of  the  tube  is 

A = F/rcv  - 8.488  X 1 06/1 .79  X 10®  = 47.42  sq  in. 


The  radius  is  given  by  Eq.  4.3  6-2 


47.42  = 13.74  in. 
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The  moment  of  inertia  of  the  section  is  given  by  Eq.  4.3.6-3. 


1 = = ^(13.74)<  = 4773  in.4 

The  radius  of  gyration  is  given  by  Eq.  4. 3. 6-4 

p = \ 1 1K  = v 4773/47.42  = 9.71  in. 


The  effective  length  of  the  column  is  given  by  Eq.  4.3.6-6 
('  = K < = (0.7)  (150)  = 105  in. 


The  critical  buckling  stress  is  given  by  Eq.  4. 3.6-5 


rCR  = 179.000-27.95  (Clp)2 
= 179.000-27.95  (105/9.71  )2 
*=  175,000  psi 

The  critical  cross  sectional  area  is  48.50  sq  in. 


The  thickness  of  the  tube  is  0.562  in. 


The  design  algorithm  for  the  vertical  stiffened  plate  structure  follows.  The  plate  is  120  x 120  in 
and  carries  4 x 106  lbs  in  compression  and  8x10*  lbs  in  shear.  The  plate  thickness  required  to 
carry  the  shear  load  is 


t-  = — — 
s 120  r, 


8 X 10s 


sy 


(120)  (73  X 103) 


= 0.913  in. 


The  plate  thickness  required  to  carry  the  compressive  load  is 


tr  = — 

c 120  t. 


cy 


4 X 106 

(120)  (179  X 103) 


= 0.186  in. 


The  actual  plate  thickness  required  to  carry  the  combined  load  is  between  0.913  inches  and 
0.913  + 0.186  = 1.099  in.  and  is  found  by  trial  and  error  by  selecting  a thickness  and  then 
checking  to  see  if  it  satisfies  the  interaction  Eq.  4.3. 6-8 


Let  t = 1.05  in. 


The  working  compressive  stress  is 


f - _ 

c 120 1 


4 X 106 
(120)  (1.05) 


= 31 ,750  psi 


The  working  shear  stress  is 


f 

s 120 1 


8 X 106 
(120)  (1.05) 


= 63,500  psi 
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The  critical  buckling  stresses  are  calculated  using  Eq.  4.3.6-10 
rCR  = K E (t/b)2 

tCCR  = (635)  <29  0 X 10*)  0-0 5/24)2  = 352,000  psi 
so  set  r = 179,000  psi 

CCH 

r = (8.00)  (29.0  x 106)  (1.05/24)2  = 444,100  psi 

SCR 

so  set  r = 73,000  psi 

SCR 


The  interaction  Eq.  4. 3.6-8  is 

*V-5  + Rc  * 1 

(VtSCR^  ' 5 + ^c^ccr)  ' ^ 

( 63,500  y 5 / 31,750  \ 

V 73,000  / ' 179,000  ' U 


.989  < 1 


; | 


Therefore,  a plate  thickness  of  1.05  in.  will  satisfactorily  carry  the  combined  load. 

A.4.2  COMPUTER  PROGRAM  LOGIC  CHARTS 


A.4.2.1  AIRCRAFT  PARAMETER  LOGIC:  The  computer  program  inputs  mission  require- 
ments and  wing  planform,  iterates  through  the  design  method  contained  in  Section  4.  and 
outputs  resultant  aircraft  parameters.  The  program  variables  are  shown  in  Figure  A.4.2. 1-1. 
The  program  logic  is  shown  in  Figure  A.4.2. 1-2.  The  following  symbols  are  used  in  the  logic 
diagram: 


ALT 

Altitude  (ft) 

Aspect  ratio 

^LMAX 

Maximum  lift  coefficient 

dTO 

Takeoff  distance  (ft) 

Gw 

Aircraft  gross  weight  (lbs) 

K 

1000 

L/D 

Lift  to  drag  ratio 

OUT 

Output 

R/C 

Rate  of  climb  (fpm) 

t/c 

Thickness  ratio 

V 

Velocity  (ft/sec) 

W/S 

Wing  loading  (Ibs/sq  ft) 

WT 

Weight  (lbs) 

A 

Taper  ratio 
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INPUT 

VARY 

OUTPUT 

CRUISE  SPEED 
CRUISE  ALTITUDE 
ASPECT  RATIO 
TAPER  RATIO 
THICKNESS  RATIO 
TAKEOFF  DISTANCE 
FUSELAGE  AREA 
MAX  LIFT  COEFFICIENT 

GROSS  WEIGHT:  800,000  - 2,400,000  LBS 
WING  LOADING:  20  - 160  LB/SO  FT 
ENGINE  NUMBER:  2 • 10 
REACTOR  TYPE:  GAS  or  LIO  METAL 
REACTOR  NUMBER:  1 or  2 
TAKEOFF:  JP-4  or  NUCLEAR 
SPEED 
ALTITUDE 

PAYLOAD 

STRUCTURE  WEIGHT 
ENGINE  WEIGHT 
CREW  ft  MISC.  WEIGHT 
TAKEOFF  THRUST 
CRUISE  POWER 
REACTOR  WEIGHT 
OPTIMUM  WING  LOADING 
WING  AREA 
UD 

Figure  A.4.2.1-1  Aircraft  Design  Parameters 


1: 
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Figure  A.4.2.1-2.  Aircraft  Parameter  Logic 
(Page  1 of  5) 
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Figure  A.4.2.1-2.  Aircraft  Parameter  Logic 
(Page  4 of  5) 
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Figure  A.4.2.1-2.  Aircraft  Parameter  Logic 
(Page  5 of  S) 

A.4.2  2 WING  ANALYSIS  LOGIC:  The  computer  program  performs  most  of  the  analysis  of 
Section  4.3.4.  It  takes  the  wing  design  criteria  and  calculates  wing  loads  and  section  proper- 
ties at  several  wing  stations.  The  program  variables  are  shown  in  Figure  A.4.2.2-1.  The 
program  logic  is  shown  in  Figure  A.4.2.2-2. 


INPUTS 

OUTPUTS 

LIFT  OF  THE  WING,  L 
WING  LOADING,  W/S 
ASPECT  RATIO.  « 

TAPER  RATIO,  A 

WING  STATION  SPACING,  Ay 

WING  WEIGHT  ESTIMATE,  W, 

LOAD  FACTOR,  LF 
ALLOWABLE  STRESS,  T 
MINIMUM  MATERIAL  GAUGE.  ( MIN 

AT  EACH  WING  STATION: 

WING  STATION,  yj 
WING  BEAM  LOADING.  W( 

SHEAR  FORCE,  Pj 

BENDING  MOMENT.  Mt 

HEIGHT  OF  BEAM  ENVELOPE,  hij 

CHORDWISE  LENGTH  OF  BEAM  ENVELOPE,  Vlj 

AREA  MOMENT  OF  INERTIA,  1, 

BEAM  CAP  THICKNESS,  tcj 
BEAM  WEB  THICKNESS,  twj 
LOCAL  CHORD,  Cj 
LOCAL  TOROUE  ARM, 

TORSIONAL  STIFFNESS.  Jj 
INFLUENCE  COEFFICIENT,  j Jdylj 

BEAM  CAP  WEIGHT,  W, 

BEAM  WEB  WEIGHT,  W2 
WING  SKIN  WEIGHT,  W3 
FLAP  WEIGHT,  W4 
TOTAL  WING  WEIGHT,  Wm 
WING  AREA,  Sw 
WING  SPAN,  byy 

TIP  CHORD,  CT  1 

ROOT  CHORD,  CR 

Figure  A.4.2.2-1  Wing  Deaign  Variables 
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A. 4.3  AIRFOIL  CHARACTERISTICS 


Selection  of  the  airfoils  was  discussed  in  Section  4.3.4.  The  NACA  653-618  airfoil  was 
selected  as  the  root  section.  The  NACA  65-410  airfoil  was  selected  as  the  tip  section.  This 
appendix  is  included  to  present  the  detailed  coordinate  definitions  of  the  airfoil  shapes  and 
the  aerodynamic  characteristics  of  the  airfoils.  Figure  A.4.3.-1  lists  the  airfoil  shape  coordi- 
nates for  both  airfoils.  Figure  A.4,3-2  shows  the  aerodynamic  characteristics  of  the  root 
section  and  Figure  A. 4.3-3  shows  the  aerodynamic  characteristics  of  the  tip  section.  For  the 
root  section,  the  slope  of  the  lift  coefficient  curve,  ao,  is  6.60  per  radian.  For  the  tip  section,  ao  is 
6.17  per  radian.  Although  ao  varies  uniformly  along  the  wing  from  the  root  to  the  tip,  the  root 
section  value  was  used  in  the  divergence  speed  calculations  because  it  yields  a slightly 
conservative  value  of  the  divergence  speed.  The  data  presented  in  this  appendix  were  taken 
from  Theory  of  Wing  Sections  (Ref.  1,  p.  434,  439,  614-615,  640-641). 


NACA  653  618 

(STATIONS  AND  ORDINATES  GIVEN  IN 
PER  CENT  OF  AIRFOIL  CHORD) 


UPPER  SURFACE 

LOWER  SURFACE 

STATION 

ORDINATE 

STATION 

ORDINATE 

0 

0 

0 

0 

0.172 

1.446 

0.828 

-1.146 

0.385 

1.776 

1.115 

-1.356 

0.839 

2.293 

1.661 

- 1.651 

2.026 

3.268 

2.974 

- 2.152 

4.462 

4.776 

5.538 

- 2.880 

6.936 

5.971 

8.064 

-3  427 

9.431 

6.978 

10.569 

- 3.876 

14.455 

8.602 

15.545 

-4.564 

19.506 

9.848 

20.494 

- 5.072 

24.574 

10.803 

25.426 

-5.433 

29.652 

11.504 

30.348 

-5.672 

34.738 

11.972 

35.262 

-5.792 

39.826 

12.210 

40.174 

- 5.784 

44.915 

12.186 

45.085 

-5.616 

50.000 

11.877 

50.000 

- 5.259 

55.077 

11.293 

54.923 

-4.723 

60.141 

10.479 

59.859 

- 4.053 

65.189 

9.482 

64.811 

-3.302 

70.219 

8.338 

69.781 

75.230 

7.075 

74.770 

- 1.705 

80.220 

5.719 

79.780 

-0.943 

85.189 

4.306 

84.811 

-0.268 

90.138 

2.863 

89.862 

0.239 

95.068 

1.433 

94.932 

0.463 

100.000 

0 

100.000 

0 

L.E.  RADIUS:  1.96 

SLOPE  OF  RADIUS  THROUGH  L.E.:  0.253 

NACA  65  410 

(STATIONS  AND  ORDINATES  GIVEN  IN 
PER  CENT  OF  AIRFOIL  CHORD) 


UPPER  SURFACE 

LOWER  SURFACE 

STATION 

ORDINATE 

STATION 

ORDINATE 

0 

0 

0 

0 

0.861 

0.628 

- 0.661 

0.607 

1.061 

0.893 

- 0.781 

1.089 

1.372 

1.411 

- 0.944 

2.318 

1.935 

2.682 

- 1.191 

4.797 

2.800 

5.203 

- 1.536 

7.289 

3.487 

7.711 

- 1.791 

9.788 

4.067 

10.212 

- 1.999 

14.798 

5.006 

19.817 

5.731 

| 

24.843 

6.290 

. -2.710 

29.872 

6.702 

30.128 

- 2.814 

34.903 

6.983 

35.097 

39.936 

7.138 

40.064 

44.968 

7.153 

45.032 

7.018 

50.000 

55.029 

6.720 

54.971 

- 2.340 

60.053 

6.288 

59.947 

-2.004 

65.073 

5.741 

64.927 

-1.621 

70.085 

5.099 

69.915 

- 1.211 

75.090 

4.372 

74.910 

- 0.792 

80.088 

3.577 

79.912 

- 0.393 

85.076 

2.729 

84.924 

- 0.037 

90.057 

1.842 

89.943 

95.029 

0.937 

94.971 

0.327 

100.000 

0 

100.000 

0 

L.E.  RADIUS:  0.687 

SLOPE  OF  RADIUS  THROUGH  L.E.:  0.168 

Figure  A.4.3-1  Airfoil  Shape  Coordinates 
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MOMENT  COEFFICIENT.  C_ 


'C/4 


SECTION  LIFT  COEFFICIENT.  C, 


Figure  A.4.3-3.  Tip  Section  Characteristics 


APPENDIX  A.6 
ENGINES 


A.6.1  RELIABILITY 

In  evaluating  the  probability  of  system  failure  where  the  systems  are  JP-4  engines, 
dedicated  nuclear  engines,  dual  mode  nuclear  engines,  or  a combination  of  JP-4  and  dedi- 
cated nuclear  engines,  the  assumptions  are: 


1)  The  probability  of  failure  of  the  engine  to  deliver  fuel  P(1)  is  s 0.1 

2)  Probability  of  engine  failure  due  to  some  reason  other  than  not  getting  the  fuel 
P(2)  s 0.1 

3)  Probability  of  the  nuclear  heat  exchanger  becoming  inoperative  for  some  reason 
P(3)  s 0.1 


These  assumptions  imply  that  the  product  of  the  probabilities  is  much  smaller  than  the 
sum.  Thus,  from  Figure  A.6. 1-1,  the  probability  of  failure  of  a JP-4  engine  is 

Pi  = P(1)  + P(2)  - P(1)P{2) 
but  P<1)P(2)  s min  |P(1),  P(2)| 


thus 


Pi  = P(1)  + P(2) 


The  figure  presents  the  system  in  terms  of  increasing  reliability  so  that  the  last  one 
presented,  i.e.,  the  combination  of  JP-4  and  dedicated  nuclear  engines,  offers  the  greatest 
reliability.  In  other  words,  this  combination  has  the  least  probability  of  total  system  failure. 


The  validity  of  the  assumption  that  P(1)  and  P(2)  are  less  than  0.1  derives  from  the 
relationship  that  the  total  engine  failure  probability  is  equal  to  P(1 ) + P(2).  Data  for  the  General 
Electric  CF6  engine  (Sections  8. 1.4.4  and  Appendix  A.8. 1.2.2)  give  a total  engine  failure 
probability  of  approximately  0.03;  therefore,  both  P(1)  and  P(2)  must  be  less  than  0.1. 

The  assumption  that  P(3)  is  less  than  0.1  is  based  on  a liquid  metal  heat  exchanger  failure 
rate  taken  from  a nuclear  systems  failure  data  handbook  (Ref.  52,  p.  21 1 ).  That  failure  rate  of 
3.2  x 10's  failure  per  hour  was  used  with  a negative  exponential  failure  probability  distribution 
with  time  equal  to  mission  time  of  330  hours  to  yield  a failure  probability  of  0.01.  Since  this 
point  calculation  was  one  order  of  magnitude  less  than  0.1 , the  assumption  that  P(3)  would  be 
less  than  0.1  was  felt  to  be  reasonable. 
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SYSTEM 


PROBABILITY  OF  FAILURE 


JP  4 ENGINE 


FUEL 

DELIVERY 

ENGINE, 
OTHER  ! 

Pill 

P<2) 

P,  ■ P(1) + P(2( 


2.  DEDICATED  NUCLEAR  ENGINE: 


1 


HEAT 

EXCHANGER 

ENGINE. 

OTHER 

PI3) 

P(2) 

P2  - P(2)  + P(3) 


3.  DUAL  MODE  ENGINE: 


BUT 


P3-P(2)  + P(1)XP(3) 
P(1)XP(3)£MIN(P(1),P<3» 


COMBINATION  OF  JP-4 


AND  DECIATED  NUCLEAR  ENGINES: 


P4-P(2)XP<1>+P(2)+P(3) 
+ P(1I  XPI3I 


BUT 


P(2)  X |P(1)+P(2)+P(3)]<P(2I 


Figure  A.6.1-1.  Reliability  of  Engine  System  Concepts  and  Combinations 
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A.6.2  DERIVATION  OF  OPTIMUM  EQUATIONS  FOR  THE  BRAYTON  CYCLE 


The  net  work  obtained  per  lb  of  working  fluid  is  the  difference  between  the  work  output  of 
the  turbine  and  the  work  required  by  the  compressor.  Assuming  constant  specific  heats: 


wn  - w,  - wc  - cp  ,,  (T„  - ru.)  - (Tu.  - T,,) 


(A.6.2-1) 


where: 


Xt  " (Pt3/Pt4> 

xc  " <fVpt,> 


k = 21_L 


Assuming  pressure  losses  are  negligible,  which  implies  that  X{  = Xc  = X,  then  the  pressure 
ratio  corresponding  to  maximum  net  work  may  be  found  by  differentiating  Eq.  A.6.2-1  with 
respect  to  X (Ref.  77,  p.  51). 

dWn  cpTt3T»t  cpTt, 
d X ~ X2  r,c  “ 0 

’»C**tTt3 


where: 


therefore: 


X = V Z tjc  ijt 

Z = Wu 

1 

rp  = l 2 VC  V\  1 2 k 


(A.6.2-2) 


An  alternative  method  for  obtaining  the  pressure  ratio  for  max  net  work  is  discussed  in 
Reference  45,  p.  194,  where  an  ideal  process  is  assumed  and  X{  = Xc  - X.  For  an  ideal 
(isentropic  process) 

evpt/  - T,/Tu  (VP,,)  * - T,/r,t 

also 

Pt2  = Pt3  and  Pt,  = Pt4 
T Tt3  Tt3  Ttl 

Substituting  this  into  Eq.  A.6.2-1 


wn=cp 
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and  assuming  T(l  to  be  constant  and  differentiating  Wn  with  respect  to  T(j,  i.e.  compressor 
outlet  temperature: 


Tt,  Tt3 


-1  = 0 


which  gives 


Tt,  = vtTt: 


(A.6.2-3) 


Once  this  temperature  isdetermined,  the  ideal  pressure  ratio  may  be  found  from  the  isentropic 
relationships  presented  on  the  previous  page. 


Besides  optimizing  a cycle  for  net  work,  a cycle  can  be  optimized  for  maximum  thermal 
efficiency,  as  shown  by  Dusinberre  (Ref.  45,  p.  31 ) and  Hosney  (Ref.  77,  p.  42).  Starting  with  Eq. 
A.6.2-1  and  noting  that  the  thermal  efficiency  of  a cycle  is  given  by 


*»th  = 


Then  for  constant  specific  heats: 


[ ,.l] 

L \/  J «l  * 


Tt,  - Tt,  [ 


Rearranging  terms: 


»tz  ( 1 -zr)  -'-(X-  D 

A TJC 


(A.6.2-4) 


To  find  ootimum  pressure  ratio  corresponding  to  maximum  efficiency,  differentiate  with 
respect  to  X and  set  equal  to  0. 

dv  du 

“1th  u dx  ~ v clX  _Q 
d X v5 


where: 


-z-(*^-+ 

V »»C 


then: 


dv  _ 1 . du  _ vl z _ ]_ 

dlT  " rfc  an°  dX  X*  i)c 
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I 


and: 


! 

i 


multiply  through  by  X2  »jc  to  obtain 

ijc  tjj  Z2  - Z X - rjj  Tjc  Z + Z - Z X2  + X2  + tjj  Z X 2 - Z X = 0 
collecting  terms 


l fo,  - 1)  Z + 1 ] X2  - 2 Tjt  Z X + [ ( 1 - tjc  ) i;t  + ijc  »jt  Z ] Z = 0 

solving  the  quadratic  for  the  pressure  ratio  which  is  defined  by 
1 


r = X ^ 
rP  * 


thus: 


'p  ijmax 


-[ 


b-\b2-4ac 
2 a 


(A.6.2-5) 


where: 


a = (»}t  - 1)  Z + 1 


b = 2q,Z 


C = l (1  - V?c)  »||  + T|c  1?t  Z 1 Z 


A. 6.3  METHODOLOGY  USED  IN  INDIRECT  CYCLE  HEAT  EXCHANGER  ENGINE  ANALYSIS 

Program  CARPET,  written  by  Capt.  R.  E.  Witherell  (Ref.  192),  a design  point  program,  was 
used  to  find  an  optimum  design  at  30,000  ft.  The  input  conditions  were: 

1)  turbine  inlet  temperature  T{<  = 1960°R 

2)  overall  pressure  ratio  Rp  = 16.8 

3)  mass  flow  of  air  through  the  engine  core  Mc  = 200  Ib/sec 

4)  pressure  loss  in  engine  A P = 6 to  16% 

The  methodology  followed  a format  as  that  shown  in  Figure  A. 6.3-1 . A simplified  diagram 
for  CARPET  is  presented  in  Figure  A.6.3-2  which  gives  an  idea  of  the  program's  format. 
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r 


INLET 


LP  COMPRESSOR  r 


HP  COMPRESSOR 


AIR  FLOWS 


BLEED  FLOWS , 

I } 

MAIN  COMBUSTOR  I 

I ! 


HP  TURBINE 


LP  TURBINE 


TAILPIPE  & A/B 


NOZZLE 


DUCT  & D/B 


PERFORMANCE  CAL. 


OUTPUT 


Figure  A.6.3-2.  Simplified  Flow  Diagram  for  CARPET 


A.6.4  SAMPLE  CALCULATION  OF  THE  HEAT  EXCHANGER  TURBOFAN  ENGINE 

The  example  problem  will  be  for  an  engine  having  the  following  specifications: 


1)  bypass  ratio  p = 4.0 
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2)  compressor  efficiency  »/c  = 0.88 

3)  fan  efficiency  = 0.88 

4)  turbine  efficiency  = 0.9 

5)  flight  speed  Mach  = 0.6,  VQ  = 596  ft/sec 

6)  altitude  = 30,000  ft 

7)  fan  pressure  ratio  Rpf  = 1.4 

8)  overall  engine  pressure  ratio  Rp  = 16.8 

9)  ram  recovery  factor  vram  - 0.9 

10)  ambient  pressure  PQ  = 627  Ibs/sq  ft 

11)  ambient  temperature  T0  = 412°R 

The  stations  are  identified  in  Figure  A. 6.4-1. 


Figure  A. 6. 4-1.  Stations  in  Nuclear  Engine 


SOLUTION: 


A.  DIFFUSER:  From  Figure  5.3  on  page  103  of  Hesse  (Ref.  69) 


= 1.148 


= 1.072 


thus,  the  conditions  at  the  face  of  the  fan  are: 


Pt2  = (1.148)  (627  lb/ft2)  = 720  lb/ft2 
Tt2  = (1.072)  (412°R)  = 442°R 


B.  FAN:  for  a fan  pressure  ratio  of  1.4 
Pt2  5 = (1-4)  (720  lb/ft*) 


P._  = 1008  lb/ft2 


The  fan  work/lb  air  is  given  by 


...  CP  Tt2  *f  (0.24)  (0.442)  (0.1004) 
f „f  0.88 


Wf  = 12.1  BTU/lb 


k 

where  Xf  = rp  - 1 and  is  in  table  form  in  (Ref.  69) 


The  fan  horsepower  is 


mfwf  (1000)  (12.1) 

HP“  0707 “0707 * 17118 


The  temperature  at  the  fan  outlet  is  given  by 


Tt2.s  ' Tt2  + 


Ju  s = 442"R  + 50.42°R 


Tt2.5  = 492  4°R 


C.  COMPRESSOR:  Because  the  overall  engine  pressure  ratio  must  be  16.8,  the  high 
pressure  compressor  (HPC)  pressure  ratio  is  12.0.  Thus 


Pt3  = <12)  pt2.5  = 12096  lb/,t2 


and  the  compressor  work  is  found  from 


wc  = 


0 24  Tt2.i  xc  _ (0.24)  (492.4)  (1.34) 


Wc  = 139.95  BTU/lb 


The  horsepovyer  required  to  drive  the  compressor  is 


HPC  - WS-  - (20°L>ll”-95>..  - 39,109 
c 0.707  0.707 


Temperature  at  the  compressor  outlet  is 


T„  - T,„  * ( 'pk  - ' 1 


T.  = 492.4  + l (12)0  295  - 1 J 

t3  0.88  ' ' 


Tt3  = 1069°R 
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D.  TURBINE:  To  find  the  turbine  exit  pressure- 

y-i 


pt5  = r *»r  po  t ( 1 + — Mo  ) ( 1 P — - ) 

,5  r o l X 2 0 ' ' »||  or  rjc  »>,  a rjf  • 

which  is  an  equation  presented  by  Hesse  (Ref.  69,  p.  272)  where 


] 


1_ 

k 


r = total  engine  pressure  ratio 
« = total  cycle  temperature  ratio  T^/T^ 


thus 


Pts  = 1522  lb/ft2 

Checking  to  see  if  the  nozzle  is  operating  supercritically: 


( 1 

1.03 

1 01 

V 1 

(0.9)  (4.43)  0.88 

(0.9)  (4.43)  0.88 

/ J 

3.5 


P 

5-5-  = 0.4133 
pt5 

which  is  less  than  the  0.495  required  for  supercritical  operation  with  a nozzle  efficiency  of 
0.98;  therefore,  the  nozzle  is  supercritical. 


Total  pressure  entering  the  turbine.  Accounting  for  the  pressure  drop  across  the  heat 
exchanger  and  the  combustor  for  a dual  mode  engine,  tho  pressure  less  will  be  assumed  to  be 
11%  (Section  6.1.5).  Then,  Pt4  may  be  found  from 

Pt4  = (1  - 0.11)  P t=  10765  lb/ft2 
so  the  pressure  ratio  across  the  turbine  is 


Turbine  exit  temperature  may  be  found  from 
Tt,  = Tt.  - ^ Tt.  t 1 - — ] 

f n 


Tts  = 1205  R 


The  turbine  work  must  equal  that  required  by  the  compressor  and  the  fan,  thus 
Wt  = Wc  + (i  Wf 

Wt  = 139.95  + 4 (12.1) 

W,  = 188.4  BTU/lb 


X- 
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The  horsepower  delivered  by  the  turbine  is 


= (200)  (188.4)  , = 53  295 

1 0.707 

Because  supercritical  flow  exits,  the  exit  pressure  is  the  critical  pressure: 


Pce  = P*  = (0.495)  (1522)  = 761  lb/ft2 


similarly  for  the  exit  temperature 


T = T*  = — - — 
ce  1.165 


= 1034°R 


E.  EXIT:  The  exit  velocity  of  the  engine  core  is  given  by 


Ve  = v'2yR  T* 


Va  = 1576  ft/sec 

The  area  exit  of  the  nozzle  is  found  from  the  law  of  continuity: 

m = p A V 

. m m R T 
pV  ~ P V 

A = (200  Ib/sec)  ( 53.3  ) (1034°R) 

(761  lb/ftz)  (1576  it/sec) 

Ace  = 9.2ft2 

or  a diameter  of 

Dce  = 3.42  ft 

Denoting  the  exit  of  the  fan  as  Station  6,  the  conditions  there  are 
pta  = Pt2^  = 1008  lb/ft2 

Tta  = Tt2.5  = 492.4°R 

The  exit  velocity  of  the  bypass  air  is  given  by  Eq.  5.33  of  Hesse  (Ref.  69,  p.  128) 
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which  gives 


V s = 862  ft/sec 


The  exit  area  of  the  fan  nozzle  may  be  found  as  follows: 


Az  = 24.14  ft2 


r2  = 3.26  ft 


Def  = 6.51  ft 


F.  THRUST:  The  net  thrust  of  a turbofan  engine  is  equal  to  the  gross  thrust  of  the  core  and 
bypass  minus  the  ram  drag: 

m/j  mc 

Fn  = <vef  ' vo)  + -f-<Vce-V 

+ (pef  ~ po>  Aef  + (pce  — po^  ^ce 

However,  Pef  = P0,  i.e.,  the  exit  pressure  of  the  bypass  air  is  equal  to  the  ambient  air  pressure 
because  it  is  expanded  to  atmospheric  pressure. 

Thus: 

Fn  = (862-596)  + (1576-596)  + (761-629)  9.2 

Fn  = 6608  + 6087  + 1214 
Fn  = 13910  lb 

A.6.5  ADDITIONAL  GRAPHS  FOR  THE  HEAT  EXCHANGER  TURBOFAN  ENGINE 

Several  graphs  for  the  heat  exchanger  turbofan  engine  with  bypass  ratios  of  6.5  and  8.0 
showing  the  affects  of  overall  engine  pressure  ratio  and  horsepower  extraction  on  the  net 
thrust  are  presented  in  Figures  A.6.5-1  through  A.6.5-4. 
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A.6.6  FLOW  DIAGRAM  OF  COMPUTER  PROGRAM  USED  FOR  THE  DIRECT  CYCLE 

A computer  program  was  written  in  order  to  study  the  direct  cycle  ducted  fan  engine 
concept.  Input  parameters  for  the  program  were: 

1)  flight  speed 

2)  efficiencies  of  the  compressor,  turbine,  and  fan 

3)  pressure  losses 

4)  altitude  conditions 

5)  temperature  and  pressures  for  the  compressor  and  turbine. 

The  program  was  designed  to  calculate  the  net  thrust  for  various: 

1)  mass  flow  of  helium  (the  working  medium  in  the  gas  turbine  generator) 

2)  fan  pressure  ratios 

Also,  it  calculated  the  thrust  for  a fully  expanded  nozzle  as  well  as  a choked  nozzle  when  the 
pressure  ratio  across  the  exit  nozzle  met  critical  conditions. 

Figure  A.6.6-1  presents  a flow  diagram  of  the  computer  program. 


Figure  A.6.6-1.  Computer  Logic  for  Direct  Cycle  Engine 
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A.6.7  SAMPLE  CALCULATION  FOR  THE  DIRECT  CYCLE 


The  following  example  will  show  the  equations  used  in  the  computer  program  discussed 
in  Appendix  A.6.6  and  how  they  were  used. 


GIVEN: 

Po  = 629  lb/ft2 

Vo  = 596  ft/sec 
Rpc  " 27  = fVpt, 


SOLUTION: 


Tt0  = 412  R 


Tta  = 2060°R 

= 0.9,  r jc  = 0.88,  r/f  = 0.85 


rhelium  ~ 1 67 
r _ i _ BTU 

cp  helium  lb-°R 

BTU 

cp  air  ~ 0 24  jj^R- 


1.  TURBINE  WORK  vper  lb) 


whore 


**t  - CP  »t  Tt,  t i-jql 

X|  - (Ptj/Pt<)  k - 0.98  <Pt!/P,,> 

W,  - (1.242)  (0.9)  (2060)  [ 1 - ] 


wt  = 742.4  BTU/lb 


2.  COMPRESSOR  WORK  (per  lb) 

7ti 

wc  = ^ — l Xc  — 1 ) 


( 1.242  j—^-)  (650’ R)  [1.49-1] 


wc  = - 


Wc  = 449.5  BTU/lb 


3.  NET  WORK  (per  lb):  This  work  may  be  used  to  drive  the  fan  or  the  circulator  needed  to 
pump  the  helium  gas  in  the  primary  and  secondary  coolant  loops. 

W„  - wt  - wc 

wn  = 742.4  - 449.5 
Wn  = 294.9  BTU/lb 

4.  FAN  WORK:  Accounting  now  for  the  mass  flow  of  helium  and  air,  the  fan  work  in  terms 
of  BTU/sec  is  given  by 

Wf  = wn  mhelium  = (295  BTU/lb)  (70  Ib/sec) 

Wf  = 20650  BTU/sec 
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1 


also 

, Vffl'i  iX  . „ 

»>f 

and  Tts  is  equal  to  442*R  (See  Appendix  6.4).  Thus  a tradeoff  between  fan  pressure  ratio  and 
mass  flow  of  air  in  the  bypass  (cuciing)  is  presented.  Selecting  the  fan  pressure  ratio  as  1.6, 
then 

m Wf*»f 

cpair  Tts  l*f  " ^ 

majr  = 1154  lb  air/sec 

Checking  the  pressure  ratio  across  the  exit  nozzle  of  the  bypass 


and 


pte  ” pte  = rpf  pts=  1155  2 


5.  EXIT  CONDITIONS:  This  engine  has  no  core  exhaust  but  rather  just  the  air  passing 
through  the  ducting.  The  temperature  and  pressure  are  calculated  based  on  expansion  nozzle 
theory  as  subcritical  flow  exits  for  this  particular  engine. 

Pef  = 629.7  psf  (static  pressure) 


T _ T = t + CP  Tts 
Ttef  Tts  'ts  + ^ 


[Xf  - 1]  = 467°R 


Pt6  " 'pf  Pts  = 11 55-2 


The  exit  velocity  may  be  found  from 


V„-  {^Te([,-(^L)k]} 

' ts 


V* . 


Vej  = 985  ft/sec 


6.  NET  THRUST : Because  it  is  an  expanded  nozzle,  the  pressure  area  term  in  the  thrust 
equation  becomes  zero,  leaving: 


F„  = 


m 


air 


thus 


F„  = 


g 

1154 

32.2 


ivef-v0] 


(389) 


Fn  = 13941  lb 
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APPENDIX  A.7 
HEAT  TRANSFER 


A.7.1  PIPING  CALCULATIONS 

The  calculations  in  this  appendix  are  b,ased  on  a single  pipe  as  shown  in  Figure  A.7.1-1 


Figure  A.7.1*d.  Concentric  Pipe  Cross  Section 


The  parameter  values  used  in  this  development  are  those  derived  in  Section  7.3.1  for  the 
supply  line  of  the  two-loop  helium  indirect  cycle.  The  mass  flow  rate  of  70  Ib/sec  establishes 
the  required  flow  area  and  radius  R,: 


P pi  wpi 


Substituting  yields: 


R,  = .385  ft 

Use  the  thiUk-walled  cylinder  equation  with  the  10,000  hr  creep  rupture  stress  of  Haynes  1 88  at 
1600°F  of  6000  psi,  and  the  secondary  loop  operating  pressure  of  1800  psi  to  define  R2. 


stress  + pressure 
stress  - pressure 


Substituting  yields' 


R2  = 0.525  ft 


This  implies  the  pipe  thickness  is: 


R2  - R,  - 0.525  - 0.385=0.140  ft  = 1.676  in. 
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Using  an  initial  assumed  thickness  of  insulation  around  the  pipe  of  1 in.,  Ra  then  equals  0 608 
ft.  Once  the  pipe  and  insulation  dimensions  are  known,  the  thermal  resistance  of  their 
combination  can  be  calculated  using 


r = = 1.309  BTU/hrft2°F  (Ref.  74,  p.  22) 

where  Kp|  and  K,  are  the  thermal  conductivity  coefficients  of  the  pipe  and  insulation,  respec- 
tively. The  next  step  is  to  define  the  thermal  resistance  of  the  heat  transfer  fluid,  helium,  on  the 
inside  of  the  pipe,  and  the  thermal  resistance  of  the  stagnant,  trapped  air  surrounding  the 
outside  of  the  pipe.  The  temperature  of  this  trapped  air  is  assumed  to  be  -70°F  for  the  reason 
explained  in  Section  7.1. 1.1.  To  calculate  these  resistances,  it  is  first  necessary  to  define  the 
absolute  viscosity  of  the  flowing  helium.  This  is  done  by 


n = P V = 3 1 X 10'5  Ib/ft-sec 


(Ref.  74,  p.  113) 


where  = absolute  viscosity  of  helium  in  the  pipe 

v = kinematic  viscosity  of  helium  in  the  pipe 

The  next  step  is  to  calculate  the  Reynolds  number  of  the  helium  in  the  pipe  by: 


Re 


2 p » Ri 


= 3.71  X 106 


(Ref.  74,  p.  115) 


This  confirms  that  the  helium  is  in  a turbulent  flow  condition.  Next  calculate  the  film  conduc- 
tance of  the  helium  in  the  pipe  using: 


h = 0.625  Kh  (Re  Pr)°  4/2  Ri 
= 59.2  BTU/hrft2°F 


(Ref.  74,  p.  180) 


The  next  step  is  to  define  the  film  conductance  coefficient  for  the  air  surrounding  the  pipe. 
This  is  defined  by  hw: 


hw  = 0.18  (Tw  - Ta)’/3 
= 0.995  BTU/hr  ft2  °F 


(Ref.  74,  p.  199) 


where 


Ta  = temperature  of  surrounding  air 

Tw  = temperature  of  outside  edge  of  insulation 


Tw  is  assumed  to  be  a maximum  of  100°F  for  this  case.  Thus  with  the  thermal  resistance  of  the 
pipe  and  the  insulation  and  the  film  conductances  of  the  inside  fluid  and  the  outside  air 
defined,  the  actual  heat  being  transferred  or  lost  radially  per  foot  of  pipe  length  can  be 
calculated  using 


Q = 


2 jt  (Ts<  Ta) 


hR, 


-*■  r + 


= 0.2  MW 


hwRs 


(Ref.  74,  p.  23) 


where  Q = the  actual  heat  lost  divided  by  the  pipe  length. 

This  then  is  one  of  the  two  parameters  that  are  used  to  define  the  insulation  thickness.  The 
criterion  is  that  if  the  value  of  Q exceeds  1/10  of  one  percent  of  the  actual  q being  transferred  to 
the  engine,  another  inch  is  added  to  the  insulation  thickness  and  the  entire  calculation  is 
reaccomplished  beginning  with  the  redefinition  of  R3. 

The  other  criterion  used  to  define  the  insulation  thickness  is  the  maximum  outside  insulation 
temperature.  To  calculate  this,  use  the  fact  that  the  temperature  changes  are  proportional  to 
the  thermal  resistances  via 


where 

ro  = 1/hw  Ra 
rT  = rQ  + r + 1/h  Ri 

This  proportion  when  solved  for  Tw  yields 

Tw  = Ta  + r-(TS4-Ta)  = 561°F 
rT 


This  value  of  Tw  is  then  compared  to  the  maximum  allowable  temperature  at  the  surface  of  the 
insulation  of  100°F  and  if  it  exceeds  that  value,  another  inch  of  insulation  is  added  and  the 
entire  calculation  is  redone. 

For  the  case  in  point  of  helium  as  the  heat  transfer  fluid,  here  is  a summary  of  the  values 
obtained  after  the  insulation  is  increased  to  the  point  where  both  criteria  are  satisfied. 

Maximum  allowable  q —0.05  mw 
Actual  radial  q lost  — 0.05  mw 
Maximum  allowable  Tw  — 100°F 
Actual  Tw  - 94° F 
Insulation  thickness  — 8 in. 


A.7-3 


APPENDIX  A.8 
SAFETY  ANALYSIS 


A.8.0  INTRODUCTION 

This  appendix  is  a supplement  to  Section  8.  Sections  A.8. 1.1  through  A.8. 1.4  deal  with 
collection  and  modeling  of  the  data  used  In  obtaining  the  failure  rates  for  Section  8.1 . Section 
A.8.2  deals  with  the  consequence  models  ofboth  the  WASH-1400  report  and  this  study. 
Section  A.8.2. 1 through  A.8. 2.4  include  the  validation  of  the  WASH-1400  model.  Section 
A.8.2.5  gives  a numerical  example  of  calculations  leading  to  one  K factor.  The  product  of  a K 
factor  and  a population  density  give  the  number  of  deaths  for  a particular  weather  condition. 

The  final  Section  A.8. 3 contains  the  computer  logic  diagrams  used  in  modeling  all  aspects 
of  Section  8. 

A. 8.1  RELEASE  PROBABILITY  DATA 
A.8.1.1  CRASH  DATA 

A.8.1.1.1  CRASH  RATES:  Aircraft  crash  rate  data  were  taken  from  Reference  183  for 
1961  through  1973.  In  all  cases,  the  data  are  a combined  figure  for  the  two  major  aircraft 
categories,  bomber  and  cargo.  The  crash  rate  data  are  shown  in  Table  A.8.1.1. 1-1  where  the 
rate  is  the  number  of  crashes  divided  by  the  number  of  flying  hours. 

The  13  estimates  of  crash  rate  shown  in  Table  A.8.1.1. 1-1  were  modeled  using 
standard  statistical  techniques  to  determine  the  underlying  probability  function.  Estimators 
for  the  mean  and  standard  deviation  were  obtained  from: 


Mean  = A = 


i N 

i-Sxi 

N M 


(A.8.1.1. 1-1) 


Standard  deviation  = S = 


L N-1 


2 (Ai  - A)* 
L=1 


1/2 


(A.8.1.1. 1-2) 


The  estimates  so  derived  are. 

mean  = 1.42  x 10‘5  crashes/hour 
std  dev  = 4.26  x 10'6 

A standard  normal  distribution  with  Z statistic 


■7  A -1.42x10-* 
4.26  x 10‘6 


(A.8. 1.1. 1-3) 


where  a = crash  rate 
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YEAR 

NUMBER 
OF  CRASHES 

NUMBER  OF 
FLYING  HOURS 

CRASH 

RATE 

1961 

77 

3,787,155 

2.03  X 

io-* 

1962 

82 

3,599,941 

2.30  X 

10* 

1963 

51 

3,511,129 

1.45  X 

10* 

1964 

58 

3,527,498 

1.64  X 

10* 

1965 

53 

3,434,601 

1.54  X 

10** 

1966 

60 

3,490,596 

1.71  X 

10* 

1967 

60 

3,530,227 

1.69  X 

10* 

1968 

59 

4,116,823 

1.43  X 

10* 

1969 

51 

3,561,022 

1.43  X 

10* 

1970 

33 

3,064,790 

1.07  X 

10* 

1971 

27 

2,768,109 

9.75  X 

10* 

1972 

33 

2,630,359 

1.25  X 

10* 

1973 

11 

1,998,098 

5.50  X 

10** 

was  hypothesized  as  the  underlying  probability  distribution  and  was  tested  using  the 
Kolmogorov-Smirnov  (K-S)  test  at  a significance  level  « = 0.20,  and  the  normal  distribution 
was  accepted.  Table  A. 8. 1.1. 1-2  shows  the  ordered  values  of  failure  rate;  their  observed 
cumulative  probability,  ^ ; the  predicted  cumulative  probability,  P;  and  the  absolute  differ- 
ence between  observed  and  predicted  probabilities,  d.  The  maximum  value  of  d is  0.124,  and 
the  critical  value  of  d from  K-S  tables  (Ref.  97,  p.  300)  is  0. 1 90.  Since  0.1 24  is  less  than  0.1 90,  the 
hypothesis  of  normality  is  accepted. 


TABLE  A.8.1.1.1-2.  K-S  TEST  OF  FAILURE  RATES 


1.45  X 10-*  0.538  0.472  0.066 


1.54  X 10  s 0.615  0.610  0.005 

1.64  X 10**  0.692  0.696  0.004 

1.69  X 10  * 

1.71  X 10* 

2.03  X 10-* 

2.30  X 10  * 


In  addition  to  the  crash  rates,  the  frequency  of  occurrence  by  location  was  deter- 
mined. Crashes  during  takeoff,  initial  climb  after  takeoff,  final  descent  for  landing,  and  landing 
were  grouped  into  a category  called  airport  crashes.  All  other  crashes  were  denoted  cruise 
crashes.  The  relative  frequency  of  these  two  categories  was  found  to  be  80  percent  for  airport 
crashes  and  20  percent  for  cruise  crashes,  or: 

P(airport  crash|  crash  occurred)  = 0.8 

P{cruise  crash|  crash  occurred)  = 0.2 

The  frequency  with  which  aircraft  crashed  with  electrical  power  failed  was  also 
determined.  Of  439  crashes,  only  three  had  electrical  failure,  or 

P(electricai  power  failed  crash)  = 6.83  x 10*3 

A distribution  of  electrical  failure  probabilities  was  formulated  by  utilizing  the  methodology 
and  error  factor  values  given  in  WASH-1400  (Ref.  185,  Appendix  III,  p.  104).  The  failure 
probability  of  6.83  x 10'3  was  taken  as  the  median  value,  and  an  upper  bound  was  found  by 
multiplying  by  an  error  factor  of  ten.  The  lower  bound  was  found  by  dividing  the  median  value 
by  an  error  factor  of  ten.  This  methodology  generates  a log-normal  distribution  of  values  and 
the  computed  upper  value  of  6.83  x 1 0 2and  the  computed  lower  value  of  6.83  x 1 0 "*were  taken 
as  the  95th  and  5th  percentiles,  respectively. 

A.8.1.1.2  CRASH  IMPACT  SPEED:  The  crash  impact  speeds  for  large  multi  engine 
aircraft  crashes  were  compiled  in  a NASA  study  (Ref.  55,  p.  10-27)  and  were  used  in  this  report 
as  a probability  distribution. 

The  crashes  were  grouped  into  categories  of  airport  crash  and  cruise  crash,  and  the 
speeds  for  both  cases  were  modeled  as  in  Section  A.8.1. 1.1.  The  impact  speed  probability 
distributions  for  both  cases  were  found  to  be  normally  distributed  with  Z statistic  as  follows: 

For  airport  crash 

7 Va-195 

L = 66 (A.8.1. 1.2-1) 


For  cruise  crash 


VA  = airport  crash  velocity 
Vc  = cruise  crash  velocity 


(A.8.1. 1.2-2) 


The  data  used  are  given  in  Table  A.8.1. 1.2-1  for 
A.8.1. 1.2-2  for  cruise  crashes.  All  impact  speeds  are  given  in 


airport  crashes  and 
units  of  ft/sec. 


in  Table 
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CRASH 

NUMBER 


TABLE  A.8.1. 1.2-1.  AIRPORT  CRASH  IMPACT  SPEEDS 


IMPACT  SPEED 

CRASH 

IMPACT  SPEED 

FT/SEC 

NUMBER 

FT/SEC 

34 

17 

186 

60 

18 

186 

132 

19 

200 

135 

20 

209 

140 

21 

220 

145 

22 

220 

160 

23 

257 

169 

24 

270 

169 

25 

270 

169 

26 

270 

169 

27 

270 

169 

28 

275 

169 

29 

290 

169 

30 

300 

170 

31 

317 

100 

TABLE  A.8.1 .1.2-2.  CRUISE  CRASH  IMPACT  SPEEDS 


CRASH 

NUMBER 


IMPACT  SPEED 
FT/SEC 


A.8.1. 1.3  CRASH  IMPACT  ANGLES:  In  addition  to  the  impact  velocity,  the  angle  of 
impact  measured  between  the  terrain  and  the  impact  velocity  vector  was  compiled  in  Refer- 
ence 55,  p.  10-27.  These  were  also  grouped  into  categories  of  airport  and  cruise  crashes.  The 
impact  angle  data  for  airport  crashes  are  given  in  Table  A.8.1. 1.3-1  and  for  cruise  crashes  in 
Table  A.8.1. 1.3-2. 
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TABLE  A.8.1.1.3-1.  AIRPORT  CRASH  IMPACT  ANGLES 


CRASH 

NUMBER 

IMPACT 

ANGLE 

CRASH 

NUMBER 

IMPACT 

ANGLE 

1 

1 

12 

8 

2 

1 

13 

8 

3 

1 

14 

10 

4 

1 

15 

5 

1 

16 

6 

1 

17 

7 

2 

18 

45 

8 

5 

19 

50 

9 

5 

20 

50 

10 

5 

21 

65 

11 

6 

— 

— 

TABLE  A.8.1. 1.3*2.  CRUISE  CRASH  IMPACT  ANGLES 


CRASH 

NUMBER 

IMPACT 

ANGLE 

CRASH 

NUMBER 

IMPACT 

ANGLE 

1 

10 

7 

60 

2 

20 

8 

80 

3 

20 

9 

80 

4 

30 

10 

90 

5 

57 

11 

90 

6 

60 

— 

— 

j These  data  were  used  in  the  SIMPAK  computer  program  and  the  distribution  statistics 

are  summarized  in  Table  A.8.1. 1.3-3. 


TABLE  A.8.1. 1.3*3.  CRASH  IMPACT  ANGLE  PROBABILITY  DISTRIBUTION 


CRASH 

ANGLE 

LOCATION 

0.05 

0.50 

0.95 

MEAN 

1 AIRPORT 

1.00 

6.00 

44.5 

14.9 

CRUISE 

18.5 

60.0 

81.0 

54.9 

I 
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A.8.1.2  AIRCRAFT  SYSTEM  FAILURE  RATES 


| 


A. 8. 1.2.1  ELECTRICAL  SYSTEMS:  Aircraft  primary  electrical  system  failure  data  are 
based  on  historical  data  from  the  C-5  data  base  and  were  provided  by  the  Lockheed  Corpora- 
tion (Ref.  1 56).  The  C-5  electrical  system  data  are  broken  into  two  broad  categories,  Generator 
and  Generator  Subsystem.  These  two  categories  make  up  one  independent  electrical  system 
and  the  failure  rates  were  given  as  8.387  x 10'5  failures/hour  for  the  Generator  and  1 .51  x 10’3 
failures/hour  for  the  Generator  subsystem.  These  rates  were  taken  as  an  upper  bound  and  the 
possibility  of  improved  reliability  was  incorporated  by  setting  a lower  bound  one  order  of 
magnitude  lower  and  skewing  a log-normal  distribution  between  the  bounds  (Ref.  185, 
Appendix  III,  p.  17). 

The  distribution  statistics  are  summarized  in  Table  A.8.1.2.1-1. 


i !; 

i 

I 

j 

i 

t 


TABLE  A.8.1.2.1-1.  AIRCRAFT  ELECTRICAL  SYSTEM  FAILURE  RATES 


COMPONENT 

FAILURE  RATE 

.05 

.50 

.95 

GENERATOR 

8.38  X 10  * 

5.89  X 10  * 

8.38  X 10  * 

GENERATOR  SUBSYSTEM 

1.51  X 10* 

1.06  X 10  * 

1.51  X 10* 

A.8.1.2.2  ENGINES:  Engine  failure  rates  were  taken  from  Reference  60,  p.  132.  These 
rates  are  current  and  projected  rates  for  the  General  Electric  CF-6  engine  and  were  given  as 
9.5  x 10'5  failures/hour  and  4 x 1 0’5  failures/hour.  These  values  were  taken  as  the  end  points  of 
a Beta  distribution  with  parameters  (4,2).  This  was  done  to  favor  the  higher  failure  rate  end  of 
the  distribution  during  Monte  Carlo  random  sampling.  The  Beta  (4,2)  distribution  is  similar  to 
the  log-normal  used  in  Reference  185,  but  is  skewed  in  the  opposite  direction.  The  distribution 
statistics  are  given  in  Table  A.8.1. 2.2-1. 


TABLE  A.8.1. 2.2-1.  ENGINE  FAILURE  RATES 


A. 8. 1.3  NUCLEAR  SYSTEM  DATA 


A.8.1.3.1  NUCLEAR  SYSTEM  FAILURE  RATES:  The  failure  rates  for  all  components 
in  the  nuclear  system  are  given  in  Table  A.8. 1.3-1.  All  failure  rates  are  distributed  with  a 
log-normal  distirbution.  Most  of  the  data  were  taken  directly  from  Reference  185. 


TABLE  A.8.1.3-1.  NUCLEAR  SYSTEM  FAILURE  RATES 


COMPONENT 
OR  FUNCTION 

FAILURE  RATE 

.05 

.50 

.05 

REFERENCE 

AUX  GEN  START 

1.00  X 10  */D 

3.00  X 10  */D 

1.00  X 10-VD 

155,  p.  19-20 

AUX  GEN  RUN 

3.00  X 10  4/HR 

3.00  X 10»/HR 

3.00  X lO-'/HR 

185,  p.  19-20 

SWITCHES  (ELEC) 

3.00  X 10"*/D 

1.00  X 10  */D 

3.00  X 10  4 

185,  p.  19-20 

MONITOR  AND  CONTROL 

1.00  X 10  '/HR 

1.00  X 10-*/HR 

1.00  X 10  '/HR 

185,  p.  19-20 

PUMP 

3.00  X 10  */HR 

3.00  X lO  '/HR 

3.00  X 10  4/HR 

185,  p.  19-20 

PIPE  (3  IN) 

3.00  X 10  ,l/H 

1.00  X 10  "/H 

3.00  X 10  */H 

185,  p.  19-20 

ACCUMULATOR 

1.00  X 10^/0 

3.00  X 10~*/D 

1.00  X 10  J/D 

185,  p.  19-20 

VALVES,  OPERATE 

1.00  X 10  */D 

3.00  X 10~*/D 

1.00  X 10  VO 

185,  p.  19-20 

VALVES,  REMAIN 

3.00  X 10  */D 

1.00  x 104/D 

3.00  x 10  4/D 

185,  p.  19-20 

MOTORS 

3.00  X 10"*/HR 

1.00  X 10  »/HR 

3.00  X 10  »/HR 

185,  p.  19-20 

EXPLOSIVE  CHARGE 

1.00  X 10  */D 

1.00X  10*/D 

1.00  X 10  4/D 

156 

CONTROL  DRUMS 

S.00  X 10  »/D 

1.00  X 10  «/D 

5.40  X 10^/D 

176 

Data  for  control  drums  were  provided  by  the  Westinghouse  Corporation  (Ref.  176), 
and  were  found  to  be  identical  with  that  given  in  Reference  52,  p.  453.  These  data  were  then 
bounded  by  the  same  error  factor  applied  to  control  rods  in  Reference  185. 


Data  for  explosive  charges  were  provided  by  Mound  Laboratories  (Ref.  156).  There 
was  no  information  at  all  upon  which  to  estimate  the  data  variance,  so  the  given  data  were 
taken  as  the  median  and  a one  order  of  magnitude  variance  either  side  of  the  median  was  used. 

Valve  data  are  from  Reference  185,  and  include  several  types  of  valves.  The  median 
value  was  checked  against  handbook  values  for  isolation  valves  in  Reference  52,  p.  354-408, 
and  was  found  to  be  identical. 

i 

i Pipe  failure  rates  are  based  on  historical  data  from  both  nuclear  and  industrial 

| applications  and  are  defined  to  be  average  values  for  rupture  of  pipe  sections  between  major 

' , discontinuities  such  as  valves,  welds,  etc.  The  averaging  process  utilizes  a constant  value  for 

pipe  failure  rate  regardless  of  pipe  length. 

A.8.1.3.2  TRANSIENT  EVENT  RATES:  The  transient  event  rate  was  taken  from  Re- 
* ference  185,  Appendix  V,  p.  56,  where  transient  is  defined  as  a condition  imposed  on  the 

reactor  coolant  system  that  results  in  a demand  for  reactor  shutdown.  The  rate  is  based  on 
historical  data  of  such  shutdowns  in  ground-based  power  generating  systems,  and  is  given  as 
log-normally  distributed  with  5th,  50th,  and  95th  percentiles  of  5.70  x 10’4, 1.14  x 10*,  and  2.28 
x 10~3,  transients  per  hour,  respectively 
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A.8.1.4  CONTAINMENT  VESSEL  IMPACT  DATA 


A. 8. 1.4.1  CONTAINMENT  VESSEL  IMPACT  RUPTURE  DATA:  Impact  tests  of  simu- 
lated containment  vessels  were  conducted  by  the  National  Aeronautics  and  Space  Administ- 
ration (NASA)  (Ref.  144)  and  the  Air  Force  Weapons  Laboratory  (AFWL)  (Ref.  75)  to  test  the 
impact  rupture  survivability.  The  tests  simulated  containment  vessels  of  approximately  20  ft 
diameter  and  approximately  500,000  lbs  impacted  against  concrete  obstacles  at  various 
speeds.  Results  of  those  tests  are  given  in  Table  A.8. 1.4-1. 


TABLE  A.8.1.4-1.  CONTAINMENT  VESSEL  IMPACT  TESTS 


TEST 

AGENCY 

IMPACT 
SPEED  FT/SEC 

RESULTS 

1 

m 

241 

NO  RUPTURE 

2 

392 

NO  RUPTURE 

3 

i 

413 

NO  RUPTURE 

4 

i 

467 

NO  RUPTURE 

5 

580 

NO  RUPTURE 

6 

641 

NO  RUPTURE 

7 

915 

RUPTURE 

8 

■ 

NO  RUPTURE 

9 

1055 

NO  RUPTURE 

10 

| 1 1 ■ 

1080 

RUPTURE 

11 

■aSaH 

1100 

RUPTURE 

An  assumption  was  made,  based  on  References  125  and  89,  that  the  failure  probabil- 
ity of  the  containment  vessel  as  a function  of  impact  speed  was  normally  distributed.  Graphi- 
cal modeling  techniques  yielded  a mean  of  980  ft/sec  and  standard  deviation  of  90  ft/sec. 
However,  a slightly  more  conservative  approach  was  used  by  setting  P(failure) = 0.33  at  V = 91 4 
ft/sec  (1/3  of  total  failures)  and  P(failure)  = 0.95  at  V = 1100  ft/sec. 


Utilizing  the  normal  Z statistic 


Zo.33  = 
Zo.95  = 


915  - S 
( r 

1100  - s 

v 


the  solution  yields  a mean  (S)  of  954  ft/sec  and  a standard  deviation  ( «■ ) of  88  ft/sec  giving  a Z 
statistic  of 

S - 954 


Z = 


88 


S = impact  rupture  speed 


(A.8.1.4-1) 
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This  distribution  gives  the  probability  of  rupture  for  a given  impact  speed.  For 
example,  if 


S = 1000  ft/sec 

then 

Z = 1.07 

and  from  normal  probability  tables 
P(rupture)  = 0.857 

In  order  to  find  the  probability  of  rupture,  the  impact  speed  must  be  known,  but  the 
impact  speed  is  a random  variable  as  shown  in  Section  A.8.1.1.2.  To  find  the  probability  of 
rupture,  the  technique  employed  in  Reference  125,  pp.  152-158,  was  used.  The  method 
computes  the  probability  of  rupture  as  a normal  distribution  with  Z statistic: 


v ‘V  + "v2 

S = mean  impact  rupture  speed 

V = mean  impact  speed 

<rs  = impact  rupture  speed  variance 

crv  = impact  speed  variance 


Utilizing  Eqs.  A.8.1. 1.2-1  and  A.8. 1.4-1  gives  the  probability  of  rupture  for  an  airport  crash: 

954  195 
V 88z  + 66* 


P(rupture|  airport  crash)  = 2x  10'12 
and  Eqs.  A.8.1. 1.2-2  and  A.8.1. 4-1  give 

, 954  - 583  1 31 

v/  882  + 269* 

P(rupture|  cruise  crash)  = 9.51  x 10'2 

It  is  to  be  noted  that  these  figures  apply  only  to  impacts  against  hard  surfaces 
(concrete).  Further  modeling  to  account  for  different  weights  and  impact  surfaces  is  per- 
formed in  Section  A.8.1. 3.6. 
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A.8. 1.4.2  CONTAINMENT  VESSEL  IMPACT  BURIAL  DATA:  Testing  was  conducted 
by  NASA  (Ref.  143)  to  determine  burial  characteristics  of  a simulated  containment  vessel  when 
impacted  against  soil.  The  results  yielded  the  following  empirical  equation: 

D = 0.0031  SN(^  ),/Z(V  - 100)  (A.8.1. 4.2-1) 

where  D = penetration  depth  (ft) 

V = impact  velocity  (ft/sec) 

W = containment  vr  ssel  weight  (lbs) 

A = frontal  area  (sq  in) 

N = nose  performance  coefficient 
S = soil  constant 

The  penetration  depth,  D,  is  measured  along  the  velocity  vector  (Figure  A.8.1. 4.2-1) 
and  the  depth  of  burial  was  defined  as: 

H = D x Sin(tf) 

0 = impact  angle 


Figure  A.8.1 .4.2-1.  Soil  Impact  Burial 

The  nose  performance  coefficient,  N,  uses  a value  of  0.7  for  a sphere  and  for  this  study 
the  soil  constant  was  chosen  to  be  7.5  (packed  soil).  For  a 20  ft  diameter  containment  vessel 
weighing  1,000,000  lbs,  the  expression  for  H is: 

H = 0.0765  x (V  - 100)  x Sin  (0)  (A.8.1. 4.2-2) 


A. 8-10 


" 


fin  which  both  impact  velocity  and  impact  angle  are  random  variables.  A depth  of  burial 
probability  distribution  was  generated  via  Monte  Carlo  simulation  using  random  selections  of 
impact  velocity  and  impact  angle.  The  results  are  summarized  in  Table  A.8.1.4.2-1. 

TABLE  A.8.1.4.2-1.  DEPTH  OF  BURIAL  PROBABILITY  DISTRIBUTION 


A.8.2  REACTOR  SAFETY  STUDY  VALIDATION 

The  WASH-1400  report  contains  a consequences  model  that  is  constructed  for  a large 
fixed  nuclear  reactor.  Since  this  model  appeared  to  be  applicable  in  most  parts  to  a movable 
reactor,  it  was  decided  to  validate  the  model.  After  each  of  the  equations  was  researched  as  to 
its  source,  one  particular  ground  reactor  accident  sequence  was  followed  through  to  deter- 
mine its  associated  number  of  immediate  deaths.  The  sequence  was  divided  into  four  basic 
sections:  first,  determine  the  amount  of  radioactivity  released:  second,  model  the  dispersion 
of  the  release;  third,  a population  model  to  estimate  the  exposed  people;  last,  to  approximate 
the  number  of  deaths  using  a medical  effects  model.  The  particular  sequence  chosen  was  the 
one  resulting  in  the  maximum  number  of  deaths. 

I 

A.8.2.1  MAGNITUDE  OF  RELEASE:  The  accident  with  the  highest  number  of  fatalities  was  a 
category  2 release  for  a pressurized  water  reactor  (PWR)  (Ref.  185,  Appendix  VI,  p.  90).  This 
accident  assumes  that  the  core  melts  and  the  released  fission  products  can  readily  escape  to 
, the  outside  atmosphere.  The  products  leak  out  at  ground  level,  but  approximately  one-fourth 

of  them  are  trapped  in  the  containment  building. 

Table  A.8.2.1 -1  shows  the  activity  of  all  the  fission  products  produced  by  a 3200  megawatt 
reactor  using  UO2  fuel.  When  the  core  melts,  however,  not  all  of  these  are  released.  Only  the 
j products  that  are  in  a gaseous  state  will  escape  from  the  molten  mass  (See  Table  A.8.2. 1-2). 


All  releases  are  grouped  by  similar  magnitudes  of  release  into  nine  categories  for  the  PWR  and 
six  for  the  boiling  water  reactors  (BWR).  Category  2 release  for  both  reactors  assumes  that  a 
near  maximum  amount  of  fission  products  are  released  at  ground  level.  Each  sequence  in  a 
category  will  release  slightly  different  amounts  because  the  circumstances  for  each  are 
different.  WASH-1400  used  a computer  code  to  figure  the  exact  amount  released  for  each 
sequence.  The  percentages  used  in  this  validation  are  slightly  higher  than  WASH-1400  report 
used  for  a category  2 release  because  the  details  of  the  computer  code  were  not  available.  The 
WASH-1400  report  assumed  that  approximately  17%  of  the  fission  products  will  plate  out  on 
the  cooler  parts  of  the  reactor,  containment,  etc. 


{ 

I 
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TABLE  A.8.2.1-1.  SET  OF  INITIAL  ACTIVITY  OF  FISSION  PRODUCTS 

USED  IN  MODEL 


NUMBER 

NAME 

SOURCE  (CURIES  X 10') 

HALF-LIFE  (DAYS) 

1 

KR-85 

0.006 

3.900E  03 

2 

KR-85* 

0.26 

1.800E-01 

3 

KR-87 

0.52 

5.300E-02 

4 

KR-88 

0.76 

1.160E-01 

5 

SR-89 

1.1 

5.060E  01 

6 

SR-90 

0.052 

1.050E  04 

7 

SR-91 

1.3 

4.000E-01 

8 

Y-90 

0.052 

2.700E  00 

9 

Y-91 

1.4 

5.900E  01 

10 

ZR-95 

1.6 

6.550E  01 

11 

ZR-97 

1.6 

7.000E-01 

12 

NB-95 

1.6 

3.500E  01 

13 

MO-99 

1.6 

2.800E  00 

14 

TC-99* 

1.4 

2.500E-01 

15 

RU-103 

1.0 

4.000E  01 

16 

RU-105 

0.58 

1.800E-01 

17 

RU-106 

0.19 

3.680E  02 

18 

RH-105 

0.58 

1.500E  00 

19 

TE-129 

0.28 

4.800E-02 

20 

TE-129* 

0.10 

3.41 0E  01 

21 

TE-131* 

0.15 

1.250E  00 

22 

TE-132 

1.2 

Jpg 

23 

1-131 

0.85 

Mill 

24 

MI-131 

0.85 

8.050E  00 

25 

1-132 

1.2 

1.000E-01 

26 

1-133 

1.7 

8.750E-01 

27 

1-134 

2.0 

3.600E-02 

28 

1-135 

1.5 

2.800E-01 . 

29 

MI-135 

1.5 

2.800E-01 

30 

XE-133 

1.7 

5.300E  00 

31 

XE-135 

0.26 

3.800E-01 

32 

CS-134 

0.017 

7.520E  02 

33 

CS-136 

0.06 

1.290E  01 

34 

CS-137 

0.058 

1.100E  04 

35 

BA-140 

1.6 

1.280E  01 

36 

LA-140 

1.6 

1.660E  00 

37 

CE-141 

1.6 

3.280E  01 

38 

CE-143 

1.5 

1.370E  00 

39 

CE-144 

1.1 

2.850E  02 

40 

PR-143 

1.5 

1.360E  01 

41 

ND-147 

0.6 

1.1 00E  01 

42 

PM-147 

0.17 

9.600F  02 

43 

PM-149 

0.4 

2.200E  00 

44 

PU-238 

0.001 

3.200E  04 

45 

PU-239 

0.0001 

8.700E  06 
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TABLE  A.8.2.1-2.  AMOUNT  OF  RADIOACTIVE  ISOTOPES  RELEASED 
WHEN  FISSIONIZED  UOz  MELTS 


FISSION  PRODUCTS 

PERCENTAGE  RELEASED 

Kr,  Xe,  I,  Br, 

90 

Cs,  Rb, 

83 

Te,  Se,  Sb, 

15.1 

Sr,  Ba, 

10 

Ru,  Mo,  Pd,  Rh,  Tc, 

3 

La,  Nd,  Y,  Ca,  Pr  \ 

0.3 

Pm,  Np,  Pu,  Zr,  f 

Each  percentage  from  Table  A.8.2.1-2  was  multiplied  by  the  activity  available  from  Table 
A.8.2.1-1.  Their  total  shows  that  13.3  x 10*  curies  are  available  immediately  after  shutdown. 
This  is  0.336  of  the  activity  of  all  45  isotopes  listed  in  Table  A.8.2.1-1.  The  release  time  was 
assumed  to  be  one  hour. 


Using  a simplified  version  the  Way-Wigner  formula: 

A = Aot"0-2  = 13.3  x 10*(3600)*°  2 = 2.59  x 10*  curies, 
are  available  at  one  hour  but.  only  3/4  of  the  total  escapes 
so.  3/4(2.59  x 10*)  = 1 .94  x 10*  curies 

is  the  total  activity  that  escapes  to  the  atmosphere. 


A.8.2.2  DISPERSION  MODEL:  As  the  cloud  of  radioactivity  is  moved  from  the  reactor  site  by 
the  wind,  its  concentration  is  changed.  The  concentration  is  reduced  by  the  cloud  expanding 
and  by  minute  particles  falling  from  the  cloud.  The  activity  of  the  fission  products  is  constantly 
decreasing  with  time,  so  the  level  of  concentration  of  activity  in  the  cloud  is  also  time  varying. 
This  model  attempts  to  determine  the  concentration  of  the  escaping  radioactive  cloud  in 
curies  per  cubic  meter.  The  concentration  problem  is  a complex,  three  dimensional,  time 
varying  one,  but  there  are  approximations  available  for  its  solution  (Ref.  160,  p.  403). 

* 

To  calculate  the  concentration  of  the  cloud  of  radioactive  particles  that  escaped,  WASH- 
1400  used  the  Pasquill  expression: 


X(x)  " 


c Au  + vTir  rrf  x <rz(x)  u exp(h2/2rrz2) 


(A.8.2.2-1) 
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where  x(x)  = cloud  concentration  at  radial  position  x (curie-sec/m3) 

Q = source  (curies) 

( = angular  sector  width  (fraction  of  circle) 
u = wind  velocity  (m/sec) 

crz(x)  = standard  deviation  for  vertical  diffusion  (m) 
h = release  height  (m) 

cA  = effective  area  of  building  (2000  m2) — used  only  when  h = 0 and  us  1 m/sec 

This  expression  was  adapted  from  the  cross-wind  intergrated  equation  by  assuming  the  cloud 
was  uniformly  distributed  across  an  arc  of  a particular  width  (Ref.  160,  p.  1 13).  To  facilitate  the 
calculations  the  cloud  movement  was  broken  into  intervals.  Since  the  interval  length  used  by 
WASH-1400  was  not  given,  x was  incremented  in  500  meter  intervals,  as  illustrated  in  Figure 


TRAVEL  DISTANCE  (ml 


Figure  A.8.2.2-1.  Vertical  Concentration  Diatribution 


A.8.2.2.1  WIND  STABILITY:  The  angle  of  the  arc  and  the  vertical  diffusion  of  the 
cloud  are  determined  by  the  stability  of  the  wind.  Pasquill  classifies  wind  into  six  categories  of 
stability.  Figure  A. 8. 2.2-1  shows  the  vertical  diffusion  for  all  the  categories  while  Table 
A 8.2.2- 1 has  the  arc  of  lateral  dispersion  specified  in  fractions  of  a circle  and  degrees.  Table 
A.8.2.2-1  also  shows  the  probability  of  each  combination  of  wind  velocity  and  stability  classifi- 
cation. These  probabilities  come  from  averaging  the  weather  of  39  reactor  sites  located  in  27 
different  states. 


TABLE  A.8.2.2-1.  METEOROLOGICAL  PARAMETERS 


NUMBER 

STABILITY 

CLASSIFICATION 

RELATIVE 

PROBABILITY 

WIND  VELOCITY 
(M/SEC) 

SECTOR 

(FRACTION  OF  CIRCLE) 

1 

A 

1.02E-02 

8.00E-01 

1.1  IE-01 

2 

A 

1.94E-02 

1.1 1 £-01 

3 

A 

‘ 1.41  E-02 

1.1  IE-01 

4 

A 

0.96E-02 

7.00E  00 

1.11  E-01 

5 

B 

6.72E-03 

8.00E-01 

8.30E-02 

6 

B 

0.98E-02 

8.30E-02 

7 

B 

8.10E-03 

K 

8.30E-02 

8 

B 

8.60E-03 

7.00E  00 

8.30E-02 

9 

C 

1.10E-02 

8.00E-01 

5.50E-02 

10 

C 

2.95E-02 

7:B 

5.50E-02 

11 

C 

3.00E-02 

5.50E-02 

12 

C 

2.15E-02 

7.00E  00 

5.50E-02 

13 

D 

3.23E-02 

8.00E-01 

4.20E-02 

14 

D 

7.00E-02 

Wi.  . i liW 

4.20E-02 

15 

D 

7.70E-02 

4.20E-02 

16 

D 

7.20E-02 

7.00E  00 

4.20E-02 

17 

E 

5.60E-02 

8.00E-01 

2.80E-02 

18 

E 

0.98E-01 

2.80E-02 

19 

E 

7.96E-02 

2.80E-02 

20 

E 

6.22E-02 

7.00E  00 

2.80E-02 

21 

F 

7.30E-02 

8.00E-01 

2.10E-02 

22 

F 

5.91  E-02 

2.10E-02 

23 

F 

3.11E-02 

2.10E-02 

24 

F 

2.12E-02 

, 1 

2.10E-02 

25 

D(RAIN) 

0.91  E-01 

4.20E-02 
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A.8. 2.2.2  DEPOSITION:  As  the  cloud  moves  from  the  site,  it  deposits  minute  radioac- 
tive particles  along  its  path.  The  depletion  of  the  cloud  was  computed  in  the  WASH-1 400  report 
at  each  interval  using 

V 8Xk 

where  frk  = 


frk  = fraction  deposited  in  interval  k. 


2k  = \/  it/2  «rz  exp  (h2/2<rz2) 


v = the  particles  falling  velocity, 

8xk  = 500  meters  (Ref.  185,  Appendix  VI,  p.  20) 

By  assuming  that  the  release  height  is  always  zero,  then 
Zk-1.25„z 


The  falling  velocities  (v)  were  assumed  to  be  0.5  cm/sec  for  iodine  and  0.2  cm/sec  for  all  other 
particles.  These  are  approximate  values  and  are  considered  to  be  conservative  (Ref.  46,  p. 
109). 


The  total  activity  released  as  calculated  in  Section  A. 8.2.1  has  23.7%  Noble  Gas 
(Xenon  and  Krypton),  65%  iodine,  and  1 1 .3%  of  all  others.  The  Noble  Gases  do  not  fall  out  like 
the  other  two  groups.  Substituting  these  percentages  into  the  previous  equation.yields: 

, 0.65(0.005)500,  0.113(0.002)500  1.39 

jf.  + _ = — rr~ 

K Zku  Zku  "zu 

To  obtain  an  accurate  measure  of  the  total  fraction  remaining  in  the  cloud  the  following 
equation  was  used: 


FMk  = exp  { - ifrj  } 

This  gives  the  amount  of  the  cloud  over  interval  k as  the  negative  exponential  of  the  sum  of  the 
previous  fractions. 

A.8.2.2.3  DISPERSION  SUMMARY:  Following  through  the  maximum  credible  acci- 
dent, the  wind  is  given  as  classification  F (See  Table  A.8.2.2-1)  with  velocity  0.8  meters/sec. 
(Ref.  185,  Appendix  VI,  p.  90). 
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Using  Eq.  A.8.2.2-1 


= 1.94  x 108 

' 0.338  + 7.87  (0.021)  500  (7.8)  0.8 


3.76  x 10s 


curie-sec 
cubic  meter 


The  fraction  of  the  cloud  deposited  is  given  by 


,rk  = 


1.39 
7.8  (0.8) 


22% 


This  leaves 


FMk  « e 22  = 80% 


of  the  fission  products  still  in  the  cloud. 

Therefore  at  the  end  of  500  meters  of  travel  the  concentration  of  the  cloud  is  3.02  x 1 05 
curies  per  cubic  meter.  The  effective  area  of  the  building  term  (cAu)  becomes  insignificant  as 
the  cloud  moves  over  100  meters  away  from  the  building.  The  only  new  inputs  in  each  interval 
were  the  crz,  x.  and  the  decayed  value  of  Q.  The  comes  from  Table  A.8.2.2-1  and  the  Q is 
computed  using  the  cumulative  elapsed  time  from  cloud  release  to  the  start  of  the  interval. 

For  example  for  1000  meters: 


and 


<rz  = 12.8  meters 


13.31  x 10» 
(4100)°2 


1.89  x 10*  curies 


A.8.2.3  POPULATION  MODEL:  This  model  is  divided  into  a population  distribution  and 
evacuation  scheme.  It  first  computes  the  number  of  people  that  could  be  exposed,  then  it 
allows  for  some  to  escape  exposure  or  to  receive  partial  exposure. 

A.8.2.3.1  POPULATION  DISTRIBUTION:  Population  data  from  66  reactor  sites  were 
used  in  the  WASH-1400  Report  to  form  a distribution.  The  area  surrounding  each  site  was 
divided  into  16  equal  sectors.  All  the  data  were  ordered  in  descending  order  by  sector  and 
averaged.  The  averages  of  all  the  sites  are  shown  in  Table  A.8.2.3-1. 

The  worst  accident,  of  course,  assumes  the  highest  population  density  which  is  the 
top  line  on  Table  A.8.2.3-1 . The  160  people  in  the  first  two  miles  were  assumed  to  be  uniformly 
spread  throughout  the  two  miles  of  the  sector.  Since  stability  classification  F covers  only  7V2°, 
only  1/3  of  the  sector  is  covered  by  the  cloud.  So  within  2 miles  only  53  people  are  exposed  and 
out  to  5 miles  only  1 266  are  exposed . At  1 0 miles  (1 6,000  meters)  the  normal  500  meter  interval 
contains  0.78  sq  mi  with  1643  people  per  sq  mile  for  a total  of  1287  people. 
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TABLE  A.8.2.3-1.  CUMULATIVE  POPULATION  (THOUSANDS)  OF  CATEGORIES 
AT  SELECTED  DISTANCES  FOR  A 22Va°  SECTOR 


CATEGORY 

DISTANCE  (MILES) 

2 

5 

10 

20 

50 

1 (99-100%) 

0.16 

3.8 

28.0 

279 

1057.9 

2 (97-99%) 

0.2 

2 

15.6 

136 

1341.5 

3 (95-97%) 

0.45 

2.1 

11.8 

92.1 

895.8 

4 (90-95%) 

0.22 

2.8 

11.5 

52.9 

279.1 

5 (80-90%) 

0.21 

2 

7.6 

26.1 

220.0 

6 (70-80%) 

0.2 

1.1 

3.6 

13.2 

151.8 

7 (60-70%) 

0.06 

0.5 

2.3 

8 

119.5 

8 (50-60%) 

0.07 

0.5 

1.6 

5.3 

65.8 

9 (40-50%) 

0.03 

0.2 

1.1 

3.6 

39.8 

10  (30-40%) 

0.01 

0.1 

0.75 

2.3 

52.5 

11  (20-30%) 

0.1 

0.09 

0.3 

1.3 

29.9 

12  (10-20%) 

0.1 

0.04 

0.2 

0.5 

19.6 

13  (0-10%) 

0 

0 

0 

0 

0 

A.8.2.3.2  EVACUATION  MODEL:  in  most  reactor  accidents  there  will  be  at  least  an 
hour  of  warning  time  before  the  majority  of  the  fission  products  are  released.  This  time  could 
be  effectively  used  to  evacuate  a large  number  of  people.  The  expression  used  by  WASH-1400 
was  that  in  each  specified  interval  i,  the  population  was  reduced  by  the  factor  Fj  where: 

Fj  = a + (1  - a)  exp  { - X (tj  - TL  + Tp  ) 
and 

a = fraction  of  population  unaffected  by  evacuation  = 0.1 
X = measure  of  evacuation  rate  = 8.3  days-’ 
tj=  time  for  cloud  transport  to  mesh  point  i = 0.26  days 
Tj=  time  between  awareness  of  impending  core  melt  and  leakage  for  acci- 
dent type  j = 1.5  hours  or  0.063  days 

TL=  time  lag  associated  w ith  interpretation  of  data  and  issuance  of  warning  to 
evacuate  = 0.02  day. 

a,  X,  and  TL  are  constant  for  all  accidents  but  Tj  depends  on  the  accident  type. 

At  the  interval  ending  at  10  miles  (16,000  meters) 

F16  = 0.1  + (0.9)  exp  (-8.3(0.26  -0.02  + 0.063)} 

= 0.17 


A.8-18 


C au  ■■ 


. 'at 


1 


The  number  of  people  in  this  interval  was  calculated  at  the  end  of  the  previous  section  as  being 
1287.  Therefore  1287  (0.17)  or  1073  people  are  left  to  evacuate  with  the  half  life  of  two  hours. 
So  two  hours  later  537  people  remain,  at  four  hours  later  268  people  remain,  etc. 


A.8.2.4  MEDICAL  EFFECTS:  This  model  was  used  to  translate  the  number  of  people  exposed 
to  varying  levels  of  radiation  into  the  statistics  of  acute  deaths,  illness,  and  property  damage. 
Acute  deaths  were  defined  as  any  deaths  in  the  30  days  following  the  accident.  WASH-1400 
considered  the  vulnerability  of  the  lungs,  Gl  tract,  and  whole  body  to  a radioactive  dose.  Of 
these,  the  whole  body  dose  was  found  to  be  the  major  contributor.  For  the  validation,  this 
study  considered  only  the  acute  deaths  resulting  from  a whole  body  dose  (WBD). 

The  WBD  includes  the  radioactive  dose  received  from  particles  inhaled  from  the  cloud, 
from  the  cloud  itself,  and  from  the  particles  deposited  on  the  ground.  The  probability  distribu- 
tion used  by  WASH-1400  to  compute  deaths  from  the  dose  received  is  shown  in  Figure 
A.8.2.3-1. 


Figure  A.8.2.3-1.  Dose  vs  Probability  of  Death 


A.6.2.4.1  INHALATION:  As  a cloud  of  minute  fission  products  passes,  everyone 
engulfed  in  the  cloud  would  inhale  an  amount  of  radioactive  material  that  depends  on  the 
concentration  of  the  cloud.  The  concentration  in  500  meter  intervals  has  been  computed,  but 
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needs  to  be  converted  to  a dose  in  REM.  Table  A.8. 2.4-1  assists  in  making  this  conversion. 
Several  values  from  this  table  were  picked  at  random  and  verified  (Ref.  46,  p.  95).  Each  value 
was  multiplied  by  the  percentage  it  contributed  to  the  total  and  then  summed.  For  example  In2 
contributes  9.6%  of  the  clouds  activity  so  its  WBD  for  a 30  day  total  would  be  1 1 REM  per  curie. 
The  total  of  all  the  isotopes  came  to  560  (REM  per  curie). 

A breathing  rate  of  2.2  x 10'4  cubic  meters  per  second  was  assumed  and  multiplied 
times  the  WBDao  total  to  complete  the  conversion  factor  of  0.12  REM  cubic  meters/curie-sec. 


A.8.2.4.2  CLOUD  SHiNE:  When  exposed  to  a radioactive  source,  which  in  this  case, 
is  a cloud,  a person  will  receive  a dose  that  is  proportional  to  the  concentration  of  the  cloud 
Using  the  values  from  "[able  A.8.2.4-2  and  the  same  procedure  as  in  Section  A.8.2.4.1,  a 
conversion  factor  for  cloud  shine  of  0.274  REM  cubic  meters/curie-sec  was  computed.  This 
was  reduced  to  one  third  its  value  due  to  the  fact  that  people  spend  the  vast  majority  of  their 
time  indoors.  The  final  cloud  shine  conversion  was  0.0913  REM  cubic  meters/curie-sec. 

A.8.2.4.3  GROUND  DEPOSITION:  The  passing  cloud  left  radioactive  particles  which 
are  assumed  to  be  evenly  distributed  over  a 500  meter  interval.  To  complete  the  dose  of  these 
particles,  the  values  from  Table  A.8.2.4-2  were  multiplied  by  their  percentage  of  the  total 
activity  and  then  summed.  This  yields  13.17  REM  square  meters/curie/hour.  This  is  cor- 
rected for  the  effect  of  shelter  by  multiplying  by  1/3  to  yield  4.39  REM  square  meter/curie-hour. 
By  multiplying  the  fraction  of  the  cloud  deposited  in  each  interval  (from  Section  A.8.2.2.2)  by 
the  curies  available  at  the  time,  one  obtains  the  initial  deposition  activity  for  that  interval.  The 
activity  per  unit  area  is  obtained  by  dividing  by  the  area  of  the  interval.  This  multiplied  by  the 
4.39  REM  square  meters/curie-hour  yields  REM/hour  for  the  ground  deposition  dose  rate. 

A.8.2.4.4  SUMMARY:  To  obtain  the  numbe;  of  deaths,  each  interval  must  be  consi- 
dered separately.  The  WBD  from  inhalation  and  clpud  shine  could  be  immediately  computed 
by  assuming  the  cloud  concentration  remained  Constant  through  the  interval  and  that  no  one 
left  while  the  cloud  was  passing.  The  WBD  from  ground  deposition  depends  almost  com- 
pletely on  the  evacuation  model.  The  people  evacuated  in  each  two  hour  period  following 
cloud  passage  were  assumed  to  receive  the  average  ground  dose  for  that  period.  This  dose 
was  then  added  to  the  WBD  for  inhalation  and  cloud  shine  for  the  total  WBD.  With  this  total 
WBD  in  REM,  Figure  A.8.2.3-1  will  give  the  probability  of  death  for  those  people  evacuated  in 
that  two  hour  period.  The  acute  deaths  for  that  two  hour  period  became  simply  the  product  of 
the  number  of  people  evacuated  and  the  probability  of  death.  Each  two  hour  period  was 
computed  until  everyone  was  evacuated;  then  all  periods  were  totaled.  This  total  is  for  one  500 
meter  interval  so  all  intervals  were  so  computed  and  totaled.  This  was  the  grand  total  of  people 
killed  by  that  particular  accident  sequence  for  those  weather  and  population  combinations. 

The  interval  ending  at  8000  meters  (approximately  5 miles)  was  chosen  to  Illustrate 
the  procedure  of  computing  the  acute  deaths  due  to  WBD. 

1 . From  the  top  line  of  Table  A.8.2.3-1 , the  people  in  the  22V2°  sector  from  2 to  5 miles  is  (3800 
- 160)  = 3640.  . 
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TABLE  A.8.2.4-1.  INHALATION  DOSE  FACTORS  (REM/CI)  X 10> 


NAME 


KR-85 

KR-85* 

KR-87 

KR-88 

SR-89 

SR-90 

SR-91 

Y-90 

Y-91 

ZR-95 

ZR-97 

NB-9S 

MO-99 

TC-99* 

RU-103 

RU-105 

RU-106 

RH-105 

TE-129 

TE-129* 

TE-131* 

TE-132 

1-131 

01-131 

1-132 

1-133 

1-134 

1-135 

01-135 

XE-133 

XE-135 

CS-134 

CS-136 

CS-137 

BA-140 

LA-140 

CE-141 

CE-143 

CE-144 

PR-143 

ND-147 

PM-147 

PM-149 

PU-238 

PU-239 


0.0 
0.0 
0.0 
0.0 
3.90E  00 
1.60E  00 
3.52E-01 
6.00E-01 


3.90E-01 
2.00E  00 
6.61  E-01 
8.40E-03 
7.71  E-01 
6.44E-02 
2.85E  00 
9.75E-02 
2.00E-02 
4.40E  00 
7.32E-01 
1.96E  00 
2.40E  00 
2.40E  00 
1.30E-01 
5.70E-01 
4.00E-02 
2.90E-01 
2.90E-01 
0.0 
0.0 

1.50E  01 
5.50E  00 
8.20E  00 
7.20E  00 
8.35E-01 
8.10E-01 
3.37E-01 
6.60E  00 


3.60E-01 
3.1  IE-01 
3.04E  02 
2.83E  02 


r 

I 


TABLE  A.8.2.4-2.  AVERAGE  DECAY  ENERGIES  AND  DOSE  FACTORS 
FOR  GROUND  DEPOSITION  AND  CLOUD  SHINE 


GROUND 

CLOUD 

NUMBER 

NAME 

(RE  MS/H  R)/(CI/m2) 

REMS/(CI-SEC/m’) 

2.  This  sector  contains  4.12  sq  mi,  so  if  the  population  is  spread  uniformly,  there  are  882 
people  per  sq  mi. 

3.  The  area  of  the  500  meter  interval  ending  at  8000  meters  is  0.196  sq  mi. 

4.  In  the  interval  of  interest  are: 

882(0.196)  = 173  people 

5.  This  number  is  reduced  thru  evacuation  by  a factor  of  Fg  = 0.36,  leaving  111  people 
exposed  to  the  cloud. 

6.  All  of  these  people  received  a W8D  from  inhalation  and  cloud  shine  of 

Conversion  Cloud 

factors  concentration 

(0.12  + 0.09)  x 988  =207  REM 

7.  The  activity  of  ihe  cloud  at  the  beginning  of  the  interval  is  obtained  using  Eq.  8.2.1 -2 


where  A0  is  3/4  ol  the  released  fission  products  and  t is  the  time  in  seconds  from  reactor 
shutdown  until  the  cloud  reaches  the  interval. 

A = 150  x 108  curies 


8.  The  fraction  of  the  initial  activity  remaining  in  the  cloud  at  the  beginning  of  the  interval  is 
34%;  however,  at  the  end  of  the  interval  it  is  32%  This  difference  gives  a fraction  dropped 
of  2%. 

9.  The  fraction  of  activity  dropped  is  therefore  1.50  x 10'(0.02)  = 3.0  x 10®  curies. 

10.  This  was  spread  uniformly  over  the  interval  by  dividing  by  the  area: 

3 0 x 10*/507,000  = 5.92  curies  per  •'q  meter. 

11.  The  initial  dose  from  the  ground  is  found  by  using  the  conversion  factor 


5.92 


curies 
sq  meter 


/ 4 gg  REM  sq  meter  \ 
* curie  hour  ' 


= 26  REM/hour 


A.8-23 


12.  The  dose  two  hours  later  is  24  REM/hour  so  the  average  dose  is: 


(26  + 24)/2  = 25  REM  per  hour 

Therefore  the  people  evacuated  in  the  first  two  hour  period  receive  a dose  from  the 
ground  of  50  REM. 

1 3.  In  the  first  two  hours,  one-half  of  the  remaining  1 1 1 people  leave.  These  56  people  receive  a 
dose  of 


207  + 50  = 257  REM. 

From  Figure  A.8.2.3-1  they  have  a 0.45  probability  of  dying  or  0.45(56)  = 25  people  of  this 
group  will  die. 

14.  In  the  next  two  hours  the  ground  dose  decreases  to  22  REM/hour  so  the  28  peoplr  that  are 
evacuated  receive  a total  dose  of  303  REM.  This  means  (.62)  (28)  = 17  people  from  this 
group  will  die. 

15.  Continuing  in  that  manner  leads  to  a total  of  63  people  who  received  a lethal  dose.  By 
stepping  through  each  interval  in  that  manner,  a total  of  2178  people  were  killed  by  the 
chosen  accident  conditions.  This  is  very  comparable  to  the  2300  deaths  in  WASH-1400,  so 
it  was  concluded  that  the  method  used  was  a valid  one. 

A.8.3  COMPUTER  LOGIC  DIAGRAMS:  Figure  A.8.3-1  shows  the  computer  logic  diagram  for 
solution  of  the  fault  trees  and  Figure  A.8.3-2  shows  the  computer  logic  diagram  for  solution  of 
the  event  trees.  Included  in  the  computer  program  was  the  SIMPAK  subroutines  for  generation 
of  probability  distributions  and  random  sampling  therefrom. 
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APPENDIX  A.9 
LIFE  CYCLE  COST 
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A. 9.1  REACTOR 

The  nuclear  reactor  production  cost  was  obtained  by  contacting  organizations  with 
reactor  expertise.  The  following  is  a copy  of  the  letter  to  these  contractors. 


This  letter  is  a follow-up  on  the  telephone  conversation  with  I had  on  20 

November  1974.  At  that  time  we  discussed  the  problem  faced  by  the  Graduate 
Systems  Engineering  class  in  costing  a proposed  lightweight  reactor  for  a nuclear 
airplane.  As  I mentioned  at  the  time,  any  help  we  receive  has  to  be  on  an  unfunded 
basis. 

Our  method  of  costing  is  to  break  the  reactor  down  Into  as  many  parts  as  feasible  at 
this  time  and  obtain  estimates  from  experts  in  the  field.  These  estimates  are  in  the 
form  of  high,  low.  and  most  probable  cost.  The  opinions  will  be  combined  statisti- 
cally to  obtain  an  estimate  of  the  system  cost. 

By  high  cost  we  mean  a cost  that  you  feel  is  realistic  and  yet  a zero  probability  of  the 
component  exceeding  that  cost.  For  the  low  cost  we  imply  a cost  that  is  realistic  but 
with  probability  of  one  of  the  price  exceeding.  The  most  likely  cost  is,  of  course, 
your  best  guess. 

We  are  aware  of  the  “risk"  involved  with  an  estimate  of  reactor  cost,  especially  with 
an  undesigned,  high  technology  reactor.  Our  approach  can  be  presented  in  our 
study  in  such  a way  that  the  individual  cost  estimate  will  not  be  associated  with  the 
expert1'  who  gave  them.  Our  study  will  list  the  estimates  that  were  provided  as  one 
portion  and  list  the  experts  as  a group. 

Since  there  is  no  haro  design  for  the  reactor,  our  description  is,  we  hope,  sufficient 
for  cost  estimates,  and  yet  does  not  constrain  the  experts  to  any  unfeasible  points. 
The  reactor  as  we  envision  it  will  be  in  the  200-400  Megawatt  (thermal)  range.  We 
hope  that  fuel  technology  will  have  progressed  to  give  us  500  w/cc  as  a core  power 
density.  The  containment  vessel  should  fit  inside  a 25  ft  diameter  aircraft  hull.  For 
safety  reasons,  we  are  leaning  toward  a helium  cooled  reactor. 
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In  closing,  I must  stress  that  this  is  an  academic  study,  not  an  official  Air  Force 
report,  and  that  there  are  no  funds  available  for  this  study  and  no  contracts  to  be  let. 
If  you  cannot  afford  to  give  us  the  breakdown,  we  would  appreciate  any  unfunded 
help  you  could  provide. 

We  do  not  envision  this  requiring  an  inordinate  amount  of  your  time,  nor  do  we 
envision  your  having  to  do  this  immediately.  Getting  the  results  prior  to  5 December 
1974  will  give  us  plenty  of  time  to  use  them. 

Attached  is  a breakdown  sheet  for  your  convenience. 


Sincerely 


RALPH  J.  LUCZAK 
Captain,  USAF 
GSE-75M,  Student 
School  of  Engineering 


Atch:  a/s 


BREAKDOWN  OF  ESTIMATES 


COMPONENT 


HIGH  MOST  PROBABLE  LOW 


shielding 

circulation  and  power  supply 

intermediate  heat  exchange 

engineering 

reactor  vessel 

ducting  and  valving 

core  reflector  and  control 
structure 

refurbishment 

allowance  for  auxiliaries 


ATTACHMENT  1 
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A.9.2  AVIONICS 


The  distribution  of  cost  for  a 0.95  reliable  system  to  accomplish  a 14  day  mission  is  given 
in  Section  9.6.  Table  A.9.2-1  illustrates  cost  distributions  for  reliable  system  other  than  0.95. 


TABLE  A.9.2-1.  AVIONICS  COST  VS  RELIABILITY 
FOR  A 14  DAY  MISSION 


RELIABILITY 

M 

ILLIONS 

K1 

0.25 

0.50 

0.75 

1.00 

0.10 

8.04 

8.73 

9.08 

9.33 

9.77 

0.30 

10.11 

11.06 

11.41 

11.75 

12.73 

0.50 

cy 

12.53 

12.98 

13.31 

13.92 

0.70 

BO 

15.04 

15.64 

15.99 

16.77 

0.90 

Kd 

18.95 

19.62 

20.14 

21.09 
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